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ABSTRACT: 

We review the effective field theories (EFTs) developed for few-nucleon systems. These EFTs 
are controlled expansions in momenta, where certain (leading-order) interactions are summed 
to all orders. At low energies, an EFT with only contact interactions allows a detailed analy- 
sis of renormalization in a non-perturbative context and uncovers novel asymptotic behavior. 
Manifestly model-independent calculations can be carried out to high orders, leading to high 
precision. At higher energies, an EFT that includes pion fields justifies and extends the tradi- 
tional framework of phenomenological potentials. The correct treatment of QCD symmetries 
ensures a connection with lattice QCD. Several tests and prospects of these EFTs are discussed. 

With permission from the Annual Review of Nuclear and Particle Science. Final version of this 
material is scheduled to appear in the Annual Review of Nuclear and Particle Science Vol. 52, 
to be published in December 2002 by Annual Reviews (http:/AnnualReviews.org). 
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1 INTRODUCTION 

1.1 Why effective theories? 

Nuclear systems have often been described as pathologically complicated. The 
forces between the constituent nucleons are strong and non-central and the rel- 
atively small binding found in nuclei results from detailed cancellations between 
much larger contributions. Adding the fact that the basic interaction among nu- 
cleons is not completely known, especially at short distances, and the problems 
involved in the numerical solution of the Schrodinger equation for systems with 
many fermions, one can understand why nuclear structure remains an unsolved 
problem after decades of intense effort. This seems even more frustrating when 
one remembers that all nuclear processes are encoded in the QCD Lagrangian 
and parameter-free predictions could, in principle, be obtained. Despite all these 
difficulties enormous progress has been made throughout the years by the use of 
models capturing different aspects of nuclear phenomena. One dissatisfying as- 
pect of these models however is their basically ad hoc nature and the presence of 
uncontrolled approximations. These models are not derived from any basic prin- 
ciple (and certainly not from QCD) and contain information coming from decades 
of trial and error hidden behind apparently arbitrary choices of some contribu- 
tions over others. Each new improvement involves the same process of educated 
guesses and one is never sure of what a reasonable error estimate would be. Ef- 
fective field theories (EFTs) are useful by providing a systematic expansion in 
a small parameter that organizes and extends previous phenomenological knowl- 
edge about nuclear processes and by providing a rigorous connection to QCD. 
They also help with more technical but important issues that have plagued nu- 
clear physics in the past, like gauge invariance, "off-shell" effects and relativistic 
corrections, by borrowing heavily from the arsenal of field theory. 

1.2 What is an effective theory? 

Most of the uncertainty in nuclear processes comes from the short-distance inter- 
actions ( ^ 1 fm) between two or more nucleons (and photons, leptons). Even 
when one is interested only in low-energy phenomena, the short-distance con- 
tributions can be important. In perturbation theory, for instance, the infiu- 
ence of short-distance physics on low-energy observables appears in the existence 
of ultraviolet-divergent integrals, that is, in the dominance of high-momentum 
modes over the small-momentum ones. Sensitivity of large-distance observables 
on short-distance physics is not an unusual situation in physics, it is in fact 
pervasive in many fields. One way of dealing with it is to model the short- 
distance physics and solve the problem from a microscopic approach. In the case 
of nuclear systems this would lead either to a calculation of nuclear processes 
directly from QCD (which is currently impossible and would be, even if possible. 
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a highly inefficient way of approaching the problem) or to the use of meson- 
exchange/quark/skyrmion/... models. Another approach is the use of effective 
theories (0)0- Before introducing the particular case used in nuclear physics, let 
us consider effective theories in general. 

Suppose we want to study the low-energy behavior of a system described by 
some theory that we will call the "fundamental" theory. In the path-integral 
formalism, we can imagine integrating over the high-momentum modes A; < A, 
where the scale A is chosen to be much larger than the momentum scale we want 
to study. The result of this partial integration over the high-momentum modes 
will be a complicated Lagrangian containing an infinite number of terms. This 
Lagrangian, called the effective Lagrangian, will also be, in general, non-local, 
but this non-locality, arising from the momenta k < A will be restricted to a 
spatial scale ^ 1/A. One can thus expand those interaction terms in a Taylor 
series on d/A, where d stands for a derivative of the fields. The coefficients of this 
expansion do not depend on the soft momenta carried by the field of the effective 
theory and describe the hard physics within the scale 1/A. They are, however, 
functions of A (the coupling constants "run"). The soft {k < A) and hard (A; > A) 
physics are factorized in the effective Lagrangian. These effective-Lagrangian 
coefficients are usually called "low-energy constants" (LECs) since they encode 
all we need to know about the fundamental theory in order to compute low-energy 
observables. Notice that, up to this point, no approximation was made and the 
effective Lagrangian contains exactly the same information as the fundamental 
one. Calculations of observables done using the effective Lagrangian will contain 
two sources of A dependence. One is the implicit dependence contained in LECs, 
the other appears in the cutoff that should be used in those computations. These 
two sources of A dependence, by construction, cancel each other. 

One may wonder what the advantage is in separating the integration over 
momentum modes in two steps. The answer to this question depends on the 
particular situation in hand. In problems where the integration over the high- 
momentum modes can be explicitly accomplished the effective Lagrangian is a 
bookkeeping device that allows us to perform approximations in a very efficient 
way. That is the case, for instance, of non-relativistic QED (§), heavy-quark ef- 
fective theory (|3|), and high-density QCD chiral perturbation theory @. In other 
cases, as the nuclear systems considered here, we will not be able to explicitly 
integrate the high-momentum modes. We can however determine the effective 
Lagrangian by a combination of self-consistency requirements and experimental 
data. We start by considering the most general Lagrangian consistent with the 
symmetries of the underlying theory. This Lagrangian contains an infinite num- 
ber of arbitrary constants. For a fixed A, different values of the LECs describe 
different underlying theories. Just one set of these values will make our low- 
energy theory reproduce the same observables as the fundamental theory. We 
then resort to an approximation scheme: we expand the low-energy observables 
in powers of the small parameter Q/A <C 1, where Q is a low-energy scale like the 
momenta of the external particles, light masses, etc. Now the factorization be- 
tween high- and low-momentum contributions comes in handy. Instead of using 

^ The term "effective theories" has also other meanings besides the specific one assumed 
in this review. "Effective Lagrangian" sometimes is meant as the one including all quantum 
corrections. Other times the term "effective theory" is used to describe any model useful at low 
energies, whether or not there is a separation of scales and a rigorous expansion in powers of 
the momentum. 
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the full effective Lagrangian with its infinite number of terms we can argue that, 
at a given order in the Q/A expansion, only a finite number of them will con- 
tribute, since the remaining terms will include many powers of 5/A ~ Q/A. This 
way we are left with a much simpler Lagrangian, with a finite (and hopefully 
small) number of coefficients that can be determined from some experimental 
data (or from the fundamental theory, if possible, or models) and used to predict 
others. Increasing the order in Q/A of a calculation will increase its precision but 
may also bring other LECs that will have to be determined by experiment. The 
precise argument connecting the order of the expansion in powers of Q/A and the 
terms in the effective Lagrangian that need to be included at that order (called 
"power counting" ) varies case by case but always include two steps. The first one 
is to estimate the size of diagrams, given the size of the LECs appearing on the 
vertices, and it is simply done by dimensional- analysis arguments. The second 
one is to estimate the size of the LECs themselves. For that we first determine 
their running, that is, their dependence on A, by requiring physical observables 
to be A independent (at the order in Q/A we are working). The information 
about the evolution of the LECs is not by itself enough to determine them since 
we do not know their initial conditions. Although for some particular value of A 
one LEC might be passing through zero, this is very unlikely. We assume that a 
typical size for a LEC C(A) is C(A) ~ C(2A) - C(A), that is, the LECs should 
have the same order of magnitude as the size of their running. In perturbative 
settings this principle amounts to little more than dimensional analysis, and is 
known as "naive dimensional analysis" (|5|). Strictly, this provides only a reason- 
able lower bound, so one should be aware of possible violations of this principle. 
This estimate, of course, is used only in arguing that some terms in the eff^ec- 
tive Lagrangian will have a negligible effect and can be dropped. The values of 
the LECs actually kept in the calculation are determined by experimental data. 
Notice that, for a given set of symmetries and low-energy degrees of freedom, 
there is no guarantee that the effective Lagrangian can be truncated at any order 
in Q/A, that is, there is no guarantee that a consistent power counting can be 
found. 

The version of the EFT method sketched above is sometimes called the "Wilso- 
nian" effective theory. Another version of the same basic idea, identical in spirit 
but differing in detail is given by the "continuum" effective theory (^). There, 
we construct the effective theory in such a way as to reproduce the same vertices 
and propagators as the full theory at low energy. The two theories differ in the 
ultraviolet region but this difference can always be absorbed in the values of the 
LECs. The technical advantage over the Wilsonian approach resides on being able 
to integrate over all momenta (used in conjunction with dimensional regulariza- 
tion), and not only over k < A, which makes it simpler to maintain gauge, chiral 
and spacetime symmetries, and to avoid power-law divergences that sometimes 
complicate power counting. 

1.3 How? 

EFTs can be used in a few different ways in nuclear physics. Historically, the 
first one was to set the separation scale A around the p-meson mass and keep 
as low-energy degrees of freedom the pions and the nucleonsQ (and maybe the A 

^ Since energies are measured from the ground state with a given baryon number, slow nucle- 
ons, despite carrying large rest-mass energy, should be considered low-energy degrees of freedom. 
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isobars), as well as photons and leptons (^,[^,§,^^,11). This approach builds on 
and extends the success of Chiral Perturbation Theory (ChPT) in the mesonic 
and one-baryon sectors. It shares with nuclear potential models the fact that 
it describes non-relativistic nucleons interacting through a potential, but it also 
brings a number of ingredients of its own, such as a small expansion parameter, 
consistency with the chiral symmetry of QCD, and systematic and rigorous ways 
of including relativistic corrections and meson-exchange currents. 

Another way of applying EFT ideas in nuclear physics is made possible by 
the existence of shallow bound states, that is, binding energies much below any 
reasonable QCD scale (|l|,|l|,||,|l5l). We can then set A around the pion mass and 
keep as low-energy degrees of freedom only the nucleons (and photons, leptons). 
At least in the case of two- and three-body systems the bound states will be within 
the range of validity of this simpler theory. This "pionless" effective theory can 
be considered as a formalization and extension of the old effective-range theory 
(ERT) (^) and the work on "model-independent results" in three-body physics 
(p^. The new features, besides the existence of a small parameter on which 
to expand, appear in a number of new short-distance contributions describing 
exchange currents and three-body forces, as well as in relativistic corrections, that 
are transparent in this approach. An extra bonus is the possibility of deriving 
analytic, high-precision expressions for many observables that previously could 
only be obtained after non-trivial numerical work. 

In Sects. |2| and ^ we will review these two approaches in few-nucleon systems, 
emphasizing qualitative aspects of recent developments. In Sect. ^ we present an 
outlook, including other approaches that are being developed for larger nuclei. 
Some reviews have already appeared covering applications of EFT ideas to nuclear 
physics, with different emphasis from the present one (18). Many developments 
of the last couple of years are described in Ref. (^) . 



2 EFT WITHOUT EXPLICIT PIONS 



2.1 The two-nudeon system and the non-trivial fixed point 
2.1.1 Two-nucleon scattering 

Let us now apply the ideas outlined in the previous section to the specific case of 
two nucleons with momentum k below the pion scale k < m-,^. More details can 
be found in Ref. (|l5|). We start by writing the most general Lagrangian involving 
only two nucleons (electroweak external currents will be included later). A system 
with two nucleons with zero angular momentum L = can exist in a spin singlet 
(^So) or spin triplet i^Si) state so there are two independent interactions with 
no derivatives, 

£ = N\id^ + Z!, + . . .)Ar _ C^t{N^PtNf - C^N^ P,Nf + . . . , (1) 

where 



PI = 
= 

^ s 



(2) 
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Figure 1: Graphs contributing to the LO A^A^ scattering amphtude. 



are the projectors in the triplet and singlet spin-isospin states (it's act on spin 
space, r's on isospin space), M is the nucleon mass and N the nucleon field. The 
dots in Eq. (|^) stand for terms with more derivatives that, as we will argue later, 
will be subdominant. 

The nucleon-nucleon (NN) scattering amplitude can be written in terms of the 
phase shift 5 as 

47r 1 



M k cot 6 -ik 
47r 1 



M -± + r^k^ + ...-ik 



(3) 



It can be shown that for potentials of range ^ R[R,^ \ /m-„: in our case), k cot 5 is 
an analytic function around A; = and that it has a cut starting at /c^ ~ 1 / , so 
it is well approximated by a power series as shown in the last line of Eq. (^. The 
parameter (?"os) is called the singlet scattering length (singlet effective range). 
For notational simplicity we specialize for now on the spin singlet channel. 

The graphs contributing to A^A^ scattering generated by the Lagrangian in 
Eq. (|^) are shown in Fig. (|^). The L-loop graph factorizes into a power, 

L— loop graph ~ (cA — ik)^ , (4) 

each one containing a linearly divergent piece and the unitarity cut ik (in the 
center-of-mass system with total energy k'^/M). The loop integral is linearly 
divergent and the coefficient c is dependent on the particular form of the regulator 
used, that is, the particular form the high-momentum modes are separated from 
the low-momentum ones. Using a sharp momentum cutoff, for instance, we have 
c = 2/ IT, using dimensional regularization (DR), c = 0. The sum of all graphs in 
Fig. (|l|) is a geometrical sum giving 

T ^ (5) 

M-J^ + cA-^k■ 

We see then that terms shown explicitly in Eq. (||) reproduce the first term of 
the effective range expansion. The addition of terms with more derivatives will 
reproduce further terms in the effective range expansion. 

We can learn some important lessons from this simple calculation. Let us 
consider two separate situations. 

Natural case: For a generic potential with range R, the effective-range pa- 
rameters typically have similar size a ~ ro ~ -R. Using DR, Co can be chosen to 
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be Co = 4:TTa/M (this choice is caUed minimal subtraction). The effective theory 
is vahd for /c < and, in this range, T can be expanded as 

.a + ika'+(^^ + a'y^ + ..^. (6) 

Since Cq ~ a, there is a one-to-one correspondence between the order in the ka 
expansion, the number of Cqs vertices and the number of loops in a graph. The 
leading order (LO) is given by one tree-level diagram, the next-to-leading order 
(NLO) by the one-loop diagram, next-to-next-to-leading order (N^LO) by the 
two-loop diagram involving Cq and one tree-level diagram with a two-derivative 
vertex (not shown in Eq. (||)), and similarly for higher orders. We have then a 
perturbative expansion, even though the microscopic potential can be arbitrarily 
strong. If one uses a cutoff regulator the situation is slightly more complicated. 
Choosing A ~ 1 /R ~ 1/a we note that the most divergent piece of the multi-loop 
graphs is as large as the tree-level graph and must be resummed to all orders, 
while the energy-dependent part containing powers of ikCo is suppressed. The 
pieces that need to be resummed at leading order merely renormalize the constant 
Cq. The one-to-one correspondence between the order in the ka expansion and 
the number of loops is lost in any but the DR with minimal subtraction renor- 
malization/regularization scheme. The technical advantages arising from the use 
of DR and renormalization theory in this perturbative setting was used in the 
study of dilute gases with short-range interactions in Ref . ( ^0|) . 

Unnatural case: In the nuclear case the scattering lengths of the two S-wave 
channels are much larger than the range of the potential. The ^5*0 (neutron- 
proton) scattering length as is = —23.714 fm and the scattering length of 
the triplet (deuteron) channel ^Si is at = 5.42 fm, corresponding to momentum 
scales of l/os = 8.3 and l/oj = 36 MeV, respectively. Those scales are much 
smaller than the pion mass, m,r ~ 140 MeV, defining the range of the nuclear 
potential^. Actually, we will see later that the potential due to pion exchange is 
too weak to describe the low-energy phase shifts, and the physics corresponding 
to the large scattering lengths occurs at the QCD scale Mqcd ~ 1 GeV, what 
makes the discrepancy between nuclear and QCD scales even more startling . 
The origin of the fine-tuned cancellations leading to the disparity between the 
underlying scale and the 5-wave scattering lengths (and deuteron binding energy) 
is presently unknown. It does not appear in any known limit of QCD like the 
chiral limit (m^ 0) or large number of colors {Nc — > oo). We will just assume 
that this cancellation happens, track the dependence of observables on the new 
soft scale l/a<j,t and perform our low-energy expansion in powers of kR <C 1 while 
keeping the full dependence on kas^t ~ 1- The singlet NN scattering amplitude, 
for instance, will be expanded as 

k'^a'irQs 1 \ 



1 -|- ikag 2 (1 -|- ikag 

It is a little challenging to reproduce an expansion of this form in the EFT. If 
one uses a momentum cutoff, for instance, the constant Cqs has to be chosen to 

^Alkali atoms used in cold atomic traps frequently have scattering lengths much larger than 
their sizes. They can be made even larger by the use of a carefully tuned external magnetic 
field (Feshbach resonances). All the ideas and formalisms developed to deal with this fact in the 
nuclear domain can and have already been used to study the physics of atomic traps (pil,p2). 
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be Cos = (4vr/M)(l/as + cA). The one- loop graph is then suppressed compared 
to the tree- level one by a factor ~ MkC^g ~ A;/A and one would naively imagine 
that the leading-order contribution is given solely by the tree- level graph. But 
there are cancellations between the graphs in Fig. (||) and all these graphs need 
to be taken into account to reproduce the expansion of r above (^^. In the 



A^A^ scattering case considered here it is not difficult to see which graphs have to 
be included at each order, but in more complex situations this can be extremely 
tricky. A more convenient way to proceed is to use a renormalization prescription 
that shifts contributions from high-momentum modes to the LECs in such a way 
as to eliminate this "accidental" cancellations between different diagrams. One 
can determine which diagrams contribute at each order on a diagram-by-diagram 
basis (manifest power counting). One way to do that is to use DR with a "power 
divergence subtraction" (PDS) (p^)^. In this scheme, we add and subtract to 
the denominator of the bubble sum in Eq. (|5|) an amount M/i/47r, where ^ is an 
arbitrary scale, and absorb the subtracted term in a redefinition of the constant 
Cos(y^), that now is a function of /x. We have for the LO amplitude 

r_ 47r 1 

The constants Cos(^) is now chosen to be 

47r 1 

CM = -TT- , (9) 

M u 

as ^ 

in order to reproduce the LO piece of the expansion in Eq. (^). The explicit 
dependence on cancels against the implicit dependence contained in Cos(//). 
The point of this rearrangement is that if ^ is chosen so that ^ ~ l/a^, Cqs{i^) ~ 
Att/M(j, and the contribution of all diagrams in the bubble sum are of the same 
order, justifying the need to resum them. Let us see how this works in some 
detail. Denoting the soft scales l/a^, /i and k collectively by Q, the tree-level 
diagram is of the order Cqs ~ 47r/M/u. The one-loop graph contains two pow- 
ers of Cos, two nucleon propagators, each one counting as \/{k'^/M) ~ M/Q"^, 
and a loop integral with three powers of momentum (~ Q^), one of energy 
(~ Q'^/M) and the usual factor of l/47r from the loop integration, for a total 
of {ATT/MQf{M/Q'^fQ^/M ~ Att/MQ. Thus the one-loop diagram has the 
same size of the tree-level graph. The same occurs for the remaining diagrams 
and they all have to be resummed. It is interesting to note that this reshuffling of 
contributions between the divergent loop and the LECs amounts to subtracting 
the poles 1/{D — 2), where D is the number of spatial dimensions, that would 
exist in two space dimensions. One can easily go to higher orders and include 
terms with derivatives in the Lagrangian. A simple calculation (again, subtract- 
ing the pole occurring in two spatial dimensions) leads to expressions for all the 
LECs in terms of the effective-range parameters (and of the arbitrary scale 
For instance, denoting by C2n the coefficient of operators with 2n derivatives, 

C,_. = l^I^if 1 V, (10) 



M 2 



* Other schemes also solve this problem (b 
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C4S = ij ^ ^ +^ ^ , (11) 



where rig is the coefficient of the third term of the effective-range expansion 
("shape parameter"). 

The /3-function describing the evolution of the dimensionless coupling cqs ^ 
-M/iCos/47r is 

d 

fJ-Q-COsif^) = COs(/U) (1 - COs(/U)) (12) 



Note the existence of two fixed points (Pt^,|25D, the trivial (perturbative) one at 
cqs = and a non-trivial one at cqs = 1. For fj, <C l/|as|, as appropriate to 
the natural case discussed above, cqs is close to the trivial fixed point. Diagrams 
involving more Cqs vertices are suppressed by powers of cqs <C 1 and the system is 
perturbative. The value cqs ~ fJ-o-s corresponds to the naive-dimensional-analysis 
one and the effects of the Cqs operator become smaller at lower energies (the 
operator is irrelevant). On the other hand, for values of ^ ~ l/|as| or larger, as 
adequate to the fine-tuned case discussed here, the flow is close to the non-trivial 
fixed point. Since cqs ~ 1, the addition of more Cqs vertices is not suppressed and 
all graphs containing only this vertex should be resummed. The dimensionless 
coupling cos goes from the naive-dimensional-analysis value cqs ~ /xa^ to cqs ~ 1 
and its effects do not go away in the infrared {marginal operator). 

Since the ^Si scattering length is also unnaturally large (and consequently the 
deuteron is unnaturally shallow), the same power counting used in the singlet 
channel applies also to the triplet channel. A^A^ scattering in this channel is 
more complicated because nuclear forces, being non-central, mix it with the ^Di 
channel. There are new operators, starting with two derivatives, describing this 
mixing. Their coefficients are determined from an expansion of mixing angle 
analog to Eq. (0). Also, the LECs are usually determined by matching to an 
effective range expansion centered around the deuteron pole, as opposed to the 
one centered around /c = as is done in the sing let channel. The ^Si NN 
amplitude is parameterized as 

47r 1 
^-M_^+pM + ..._,fc ^^^^ 

where is the deuteron binding energy and p the effective range parameter. 

Explicit expressions for the leading terms in the Lagrangian and numerical values 
for the LECs can be found in Ref. (|l5|). 

The inclusion of external currents (photons, neutrinos, ...) is straightforward. 
All terms involving nucleons and the new fields or currents should be included, 
as long as they satisfy the symmetries of the underlying theory. In the case of 
photons, some of these terms are just those required by gauge invariance and are 
determined by minimally coupling the photon to the nucleon Lagrangian. Their 
coefficients are thus fixed by A^A^ scattering data and gauge invariance. There 
are also terms that are gauge invariant by themselves and whose coefficients 
are not determined by A^A^ scattering data alone. They represent the physics 
of exchange currents, quark effects, etc., and need to be determined through 
some extra piece of experimental data. To perform the low-energy expansion 
though, it is necessary to have an a priori estimate of their size. This estimate is 
obtained by using the fact that observables should be independent of the cutoff 
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(or fi if using DR), at least up to the order one is computing. Consider some 
two-nucleon operator of the form X = C^N^ N^Txd'^^NN , where Tx is some 
tensor in spin-isospin space. Its matrix element on two-nucleon states is given by 
the diagrams involving the operator X "sandwiched" between two two-nucleon 
scattering amplitudes and by one-loop one-body diagrams that do not involve X. 
Typically the one-body diagram is not divergent and does not introduce any /i 
dependence^ so the remaining graphs have to be /j, independent by themselves. We 
have to make a distinction now between the cases where the operator X connects 
two 5-wave states, two non-S-wave states, or one S'-wave and one non-S-wave 
state. In the first case renormalization-group invariance of the two-nucleon matrix 
element of X implies 

9 -,A-, 



where T is the LO A^A^ scattering matrix, which is independent. From that 
it follows that C^^{^) scales as ~ (^u — l/a)~^. Similarly, for the case where X 
connects one S'-wave or no S'-wave states C^{ij) scales as ~ {fj. — l/a)~^ and 
~ (/i — respectively. Using dimensional analysis to fix the powers of A we 

then have 

1 1 

^^^^^^ ~ M(l/a-;u)"A2-+i-"' ^^^^ 

where a is the number of S'-wave states the operator X can connect (either 0, 1 
or 2). 

In a nutshell, the power counting rules valid for the two-nucleon system are 
(|ll,|l|, 11,11): 



fermion line — > M/Q^ 



loop 



AirM 
d ^ Q 

do ^ qVm 



2n 



where is the coefficient of the two-nucleon interaction with 2n derivatives, 
is the coefficient of a two-nucleon operator with external current X and 2n 
derivatives, and A is the high-energy scale A ~ rriT^. 

Using this rule we can determine the contributions to A^A^ scattering at any 
given order. At LO, for instance, we have the series of diagrams shown in Fig. (^, 
with all the vertices containing no derivative. That is the only non-perturbative 
resummation necessary. At NLO we have the insertion of one C2 operator in a 
chain of Co operators. At N^LO we have two insertions of C2 and one insertion of 
C4, and so on. The resulting ^S'l phase shift, for example, is shown in Fig. @), and 
compared to the Nijmegen phase-shift analysis (PSA) (|2^ . Analytic expressions 
for the phase shifts can be found in Ref. (|l5|) . They suggest convergence for 
momenta k ^ 100 MeV, as it is reasonable for an EFT without explicit pions. 

^One exception is the two-nucleon, no-external-current C4 operator whose renormalization 
is driven by C2. This explains the apparent discrepancy between Eq. (^o|) and Eq. (p^. 
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Figure 2: ^Si NN phase shift (in degrees) as function of the center-of-mass 
momentum. The LO result is the dashed (purple) line, the N^LO the dotted 
(red) line and N^LO the thick (blue) solid curve. The dot-dashed line is the 
Nijmegen PSA. Prom Ref. (|l5|), courtesy of M. Savage. 



Electromagnetic effects in proton-proton scattering were considered in the EFT 
approach in Ref. (p7|). 

Up to this point we have considered only NN scattering, where the predictive 
power of the pionless EFT is very small. We were able however to determine 
many LECs using scattering data and understand the effects of the fine-tuning 
on the 5-wave channels. We now apply the formalism developed above to the 
computation of form factors and processes involving external currents. We will 
omit the diagrams needed to be computed and the explicit analytic expressions 
that are always available in the two-nucleon sector. They can be found in the 
literature cited. 



2.1.2 Electromagnetic form factors of the deuteron 



The matrix element of the electromagnetic current on the deuteron has the non- 
relativistic parameterization 



2Md 



1 „ . / 1 



E + E' 
2Md 



+^]^^s(9 ) [qiqj - 3? Sij ){p + p 



(17) 



where |p, e*) is the deuteron state with momentum p and polarization e*, is 
the deuteron mass, q = p' — p and the form factors are normalized such that 
Fc{0) = 1 (deuteron charge), eFM(0)/2M^ = fiD (deuteron magnetic moment) 
and Fq{0)/M^ = fig (deuteron quadrupole moment). 

At LO and NLO the computation of Fc{q^) involves only the constants Cot 
and C2t and is identical to the ERT calculation. At N^LO a one-body term 
describing the nucleon charge mean square radius {{r'^)N) appears, which is the 
first deviation from ERT (p^). Formally there are also relativistic corrections, but 
they are suppressed by powers of Q/M as opposed to Q/m-jr, and are numerically 
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small. Still they can be readily computed in EFT. Defining the deuteron charge 
mean square radius by (r^)^ = d{dFc / dq'^) , one finds 

to be compared with the experimental value {r'^)d = 4.538 fm^. 

The magnetic form factor FmIq'^) at LO and NLO is simply the electric form 
factor Fc{q^) multiplied by the isoscalar nucleon magnetic moment k„ + Kp, 
except for a new two-body term appearing at NLO without an ERT analog, 

C = -eL2ie'^^{NPiN)\NPjN)Bk + h.c. (19) 

L2 can be determined through the experimental value of the deuteron magnetic 
moment and, using this value, the momentum dependence of -Fm(9^) is then 
predicted. 

The Fq^q^) form factor involves a transition between the S- and D-wave com- 
ponents of the deuteron. At LO its value is determined by a S- to D-wave 
transition operator whose coefficient is extracted from the asymptotic D/S ratio 
of the deuteron, rjD/s- NLO there is a new two-body term 

£ = -eCQ{NPiN)^NPjN) (v'V^ - ^^^'^^ (20) 

whose coefficient Cq can be fitted to the experimental deuteron quadrupole mo- 
ment. At N^LO the only contribution comes from the finite size of the nucleon 
charge distribution {r'^)N- The value of Fq{0) is then a fit, but the momentum 
dependence is an EFT prediction. The presence of a counterterm not determined 
by A^A^ scattering at NLO indicates that fig is sensitive at the ~ 10% level to 
short-distance physics not determined by A^A^ scattering. That is probably the 
reason different potential- model calculations underpredict fiQ by ~ 5%. 

2.1.3 Deuteron polarizabilities and Compton scattering 

The deuteron is a very loose bound state so it is no surprise that its electric 
polar izability depends mostly on the large-distance part of its wave-function. 
Consequently, a model-independent prediction can be made with a high degree 
of accuracy. The electric (magnetic) scalar deuteron polarizabilities ueo (^mo) 
are defined by 

Cv = E'^ + 27raMod'U'B^, (21) 

where is a canonically normalized deuteron field. 
The EFT result for ueo up to N^LO is (|l|) 

aM 1 f 272 M73 \ . , ^ 

= (' + ^ + 1^"" + ■■■)= =^ "-""^ • 

where Dp = —1.51 fm^ is a combination of constants describing P-wave two- 
nucleon interactions and in the last line we used the fact that higher-order 
computations always build the deuteron wave-function renormalization factor 
Z = 1/(1 — jp) 0. The first two-nucleon current operator appears only at N^LO. 

®A modified method counting the effective range parameters as ros,p ~ 1/Q formalizes this 
observation (ESI). 
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Since all the inputs were determined from A^A^ scattering, we expect that any 
model describing the phase shifts well should make predictions within the esti- 
mated error. Indeed a number of models predict oeo = 0.6328 it 0.0017 fm^. 

Magnetic, vector and tensor polarizabilities, and their momentum-dependent 
analogs were not yet analyzed in the pionless EFT (see below for results using 
the pionful EFT). However, these polarizabilities determine the amplitude for 
low-energy Compton scattering off the deuteron, that was analyzed in Ref. (pOl). 
The nucleon polarizabilities appear already at NLO and Compton scattering can 
be used to extract the hard-to-measure neutron polarizabilities. 



2.1.4 Radiative capture of neutrons by protons 

The n + p — > d + 7 reaction at low energies is a key ingredient in Big-Bang 
nucleosynthesis calculations. The amplitude for this process can be expanded in 
multipoles as 

A = ieXMive'^''e*'k^e'' nP^p + sXeiv nT2T3a2a ■ tpP'e*' 



V2 



nP'[k't ■ t - t ■ ke*']p 



^eXE2S_ ^pi^f^i^ . ^ + ^ . fe« _ l^*ik . e-*]p + . . . , (23) 
v 2 3 

where n and p are the neutron and proton Pauli spinors, k is the photon momen- 
tum and e (e) is the polarization of the deuteron (photon). At low energies the 
form factor Xmiv dominates the cross section by a few orders of magnitude. Its 
computation at LO is the same as the ERT one and underpredicts the experimen- 
tal value for thermal neutron capture by 10%. This discrepancy was explained 



long ago as due to a pion-exchange current contribution (po|). In EFT (15) the 
same effect is encapsulated in the two-nucleon current 

£ = eLMiv{NP'N)\NP^N)6A3B' + h.c, (24) 

whose coefficient Lmiv can determined by the cold-capture cross section. The 
momentum dependence is then predicted. However for photon energies larger 
than a few MeV the Xeiv form factor dominates the cross section. Up to N^LO 
the computation of Xeiv involves only Cot, C2t, the P-wave interaction com- 
bination Dp encountered in the polarizability calculation, and Lmiv fitted at 
threshold. At N^LO some relativistic effects and a new term appears (31), 



C = eLEiv{NP'N)\N{ViP,A - jPiA)N)6A3E^ , (25) 

where PiA = o'2criT2TA/V^, which is fitted using data for the inverse reaction 
d+7 — > n+p. The resulting cross section for the photodissociation of the deuteron 
is shown in Fig. (^), and compared to data (^). The estimated error is of the 
order {Q/m.„f ~ 1%. These precise, analytical computations are particularly 
useful for Big-Bang nucleosynthesis codes. 

A set of polarization observables in the n + p — > 7 + d reaction was analysed 
in Ref. (pi). 



2.1.5 Neutrino-deuteron scattering and proton-proton fusion 

A complete set of reactions involving (anti-)neutrino breakup of the deuteron 
were computed in the pionless EFT approach (34). Proton-proton fusion and the 
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Figure 3: Cross section for 'j + d^n + p as function of the photon energy at 
N^LO, compared to data. From Ref. (31), courtesy of G. Rupak. 



interplay of Coulomb interactions were analysed in Ref. (^5[) . These reactions are 
essential for the understanding of solar neutrino physics since they are relevant 
for both the production {p + p — > (i + e"*" + i^e) and detection in heavy water 
detectors (through the neutral current i^ + d— > v + n+p and the charged current 
I' + d^e'^ + n + n reactions) . 

The weak interactions are described by the familiar neutral and charged current 
pieces 

-^(/z^J^ + Z+^J^+h.c), 



(26) 



with the leptonic currents Iz/j, = z^(l — 75)7^2^, /^ 



(1 — 75)7^6 and hadronic 



currents = V!^ - and = -2 sin^ 0y^,VQ^ + {I - 2 sin^ OwW^^ 
The isosinglet vector (axial) current Vq* (Aq) and isotriplet vector (axial) current 
(A^), written in terms of the nucleon fields, have contributions in the form 
of one- and two-nucleon operators. The one-nucleon operators are determined by 
the axial coupling constant = 1.26, the neutron and proton magnetic moments 
Kn and Kp, the strange contribution to the proton spin (57^755) and the strange 
magnetic moment of the proton /u^ . The two-nucleon currents contributing up to 
N^LO are 



A"- 
A' 



Lia{NP'N)^{NPaN), 
-2iL2AUok{NP'N)\NP^N)+h.c., 



Vi = 2iL,eijk{NP^N)\V^ + 'v')iNPaN), 
V'o = 2iL2{NPm)\Vj + ^j){NP'N). 



(27) 



Li = Lmiv was determined through the neutron-proton cold-capture cross sec- 
tion and L2 was determined by the value of the deuteron magnetic moment. The 
parameters L2A, {s'y^'j5s) and ^s, which are not well determined experimentally, 
have a negligible impact on the cross section (< 1%) due to the almost orthog- 
onality between initial and final states in the triplet channel. The only relevant 
unknown in a N^LO calculation is then the value of Lia- Using the estimate in 
Eq. ( p^ ) we find Lia ~ 47r/M//^ ~ 5 fm^. It is found that two potential- model 
results, one with and another without exchange-current terms, and differing by 
about 5%, can be reproduced by varying the value of Lia within this range. 
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This shows that the difference between these calculations comes from different 
assumptions about the short-distance physics. To fix this indeterminacy and, 
consequently, have predictions for the I'd reactions at the percent level the value 
of LiA needs to be determined experimentally. One possibility is the measure- 
ment of one of these reactions at one energy. The other is to extract Lia through 
another process sensitive to this term, like tritium /?-decay or muon capture on 
the deuteron. The challenges involved in this extraction from the well-measured 
value of the tritium lifetime will be discussed below. 



2.2 The three-body system and the limit cycle 

The study of three-nucleon systems using EFT is still in its infancy compared 
to its mature status in the two-nucleon sector. Only now calculations accurate 
enough are appearing for the triton-^He channel that open the door for precision 
calculations of processes involving external currents. Those processes may turn 
out to be a very important way of fixing the value of two-body LECs that are hard 
to measure in the deuteron. The new ingredients in going from two- to three- 
body systems are three-body interactions. In the absence of fine-tuning their 
typical size is determined by dimensional analysis. Since they subsume physics 
contained within the range 1/A ~ I/tott, a three-body force with 2n derivatives 
would have the typical size 

= D,J'-{N^Nf D,^ ~ (28) 

Just as it happens in the two-body force, the fine-tuning in the two-body S- 
wave channels introduces a new scale 7 ~ l/a^ that invalidates the estimate in 
Eq. (|28|). We will resort to the same argument used before to estimate the size 
of these contributions. We will demand observables to be cutoff independent 
order by order in the low-energy expansion, which determines the running of 
the three-body forces with the assumption that their typical size is set by the 
size of their running, Z)2n(A) ~ D2ri(2A) — D2n{^)- As mentioned before, it is 
unlikely that 152™ (A) is much smaller than this estimate for a particular value of 
the regulator A, or that it contains a large A-independent piece. Another way of 
looking at this way of estimating the LECs is to remember that if the inclusion 
of a particular LEG is necessary in order to have cutoff-independent results (at a 
particular order in the expansion), this LEG needs to be large enough to appear 
at that same order of the expansion. 

The doublet S-wave channel (where the triton and ^He are) has a very different 
behavior from the other channels. The physical reason is that this is the only 
channel where all three nucleons can occupy the same point in space (two spin 
and two isospin states allow for a maximum of four nucleons in the same state) . 
A system of three bosons also displays this property and is qualitatively similar to 
three nucleons in the doublet 5-wave channel. In the remaining channels either 
the angular-momentum barrier or the Pauli exclusion principle forbids the three 
nucleons to touch. One would then expect that the doublet 5-wave channel (and 
systems of three bosons) are much more sensitive to short-distance physics than 
the other channels, an expectation that we will see confirmed by further analysis. 

To avoid unnecessary complications we will present explicit expressions only 
for the case of the /S-wave three-boson system. The formulae for the nucleon 
cases in the different channels can be deduced in an analogous way and can be 
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Figure 4: LO graphs contributing to the dressed propagator of the dimeron (a) 
and to the particle/dimeron amphtude (b). 



found, for instance, in Refs. (^,38,|37). A convenient first step is to rewrite the 



action in terms of an auxihary "dimeron" field d ( |39| ) , 
Cbosons = N^ido + ^ + ...)N-Co{NNyNN-Do{NNN)^NNN + ... 

+y{dmN + dN^N^) + gdUN^N + ..., (29) 

where Cq = y^/A and Dq = —gy^jl^. The equivalence between the two La- 
grangians above can be shown by simply performing the Gaussian integral over 
the auxiliary field d. Since d has the quantum numbers of a two-particle bound 
state, it can be used as an interpolating field for the bound state. Note that the 
normalization of the field d (and the value of A) is arbitrary and all the observ- 
ables should be a function only of the combinations y^/A and g'y^/A^, not of y, 
g and A separately. 

The propagator A(p) for the dimeron field of momentum p is given by the 
sum of bubble graphs shown at the top of Fig. (^). The power counting leading 
to the necessity of summing all these graphs when the scattering length is large 
is identical to the one discussed in connection to A^A^ scattering. Taking the 
arbitrary constant A to scale as A ~ 1, we have ~ Co ~ ^TKjMQ and thus 
A(p) scales as A(|j) ~ (47r/My^)(l/(5) ~ 1. Summing all graphs, we can write 

A(p) = .\ (30) 

Let us now consider the graphs describing the scattering of one particle off the 
bound state of other two, shown at the bottom of Fig. (^). We can determine the 
impact of each graph by power counting. For each additional loop we have two 
particle propagators, one d propagator, two powers of y and one loop integra- 
tion, for a total of {M IQ'^f{\){^T: IMQ){Q^ j \t^M) ~ 1. There is no suppression 
for additional loops, and all diagrams involving an arbitrary number of particle 
exchanges contribute at LO. Graphs including the two-body operator are sup- 
pressed and start appearing at NLO. Graphs including the three-body force may 
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or may not be leading, depending of the size of the three-body forces. Since at this 
point we do not know how large they typicahy are, we wih provisionally include 
them. The chain of diagrams to be summed, contrary to the two-particle case, 
does not form a simple geometrical series and cannot be summed analytically. 
We can however consider the second line in Fig. as an equation determining 
this sum. For the bosonic case, A = 1, we have 



tkip) = K{p, k) + — + ^j^ dqq'A{q) {k{p, q) + — j t^iql (31) 

with 



K{p,q) 




p^ +pq + q'^ - ME 



where tk{p) is the scattering amplitude with all but the outgoing single-particle 
line on-shell ("half-off-shell amplitude"), k is the incoming momentum in the 
center-of-mass system, p is the outgoing momentum, and E = 3A;^/4m — 7^ is the 
total energy. The on-shell point is at j? = A;. In the case of nucleons in the quartet 
channel the same equation is obtained, but with A = —1/2 and some additional 
momentum dependence in the three-body force (since the simple momentum- 
independent three-body force does not contribute to this channel). Because all 
spins are parallel in the quartet channel, only triplet two-body interactions occur 
and the value of 7 is determined by the deuteron pole. The auxiliary field d 
carries the quantum numbers of the deuteron. In the doublet case, singlet and 



triplet two-body interactions contribute. The analog of Eq. (31) is a pair of 
coupled integral equations that, in the S'C/(4) limit[[] where the singlet and triplet 
scattering lengths are equal (or in the ultraviolet where l/a^, 7 can be discarded) 
decouple into a pair of equations like Eq. (|3l|), one with A = 1, another with 
A = —1/2. Two auxiliary fields appear, one with the quantum numbers of the 
deuteron, another with the ^5o quantum numbers. In all spin channels, equations 
for the higher partial waves are obtained by substituting the logarithm in the 
kernel by a Legendre function Qi{pq/{p^ + q^ — ME)). Eq. ( |3lD is the version 
of the Faddeev equation (see Ref. ( |4l| ) for a nice introduction) appropriate for 
contact interactions. It was first derived, by very different methods, in Ref. (Q). 
It is only for separable potentials, like the contact interactions considered here, 
that the Faddeev equation reduces to a one-dimensional integral equation. This 
simplification reduces the numerical work involved by many orders of magnitude. 



2.2.1 Most channels are easy 

Let us first consider the channels in the second group, that is, all channels but 
the doublet 5 wave. The diagrams summed by Eq. (^) are all ultraviolet finite. 
That would suggest that there is no need to include three-body forces to absorb 
the cutoff dependence since this dependence is a subleading 1/A effect. The 
(numerical) solution confirms this. Even in the absence of a three-body force the 
phase shifts are only weakly dependent on A and tk{k) has a finite limit when A — > 
00. Higher-order corrections can be included either perturbatively (as was done 

■'See Ref. (CT) for a discussion of the 5f/(4) limit in the two-nucleon sector. 
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Figure 5: k cot 5 for nd scattering in the quartet S'-wave channel as function of the 
energy. The dashed hne is the LO and the full line the N^LO result (|3^). Points 
are from a PSA ( ^3|) and a near-threshold measurement (p^). Figure courtesy of 
H.-W. Hammer. 



in the two-body sector) or non-perturbatively through the denominator in the 
dimeron propagator. This last resummation amounts to including some (but not 
all!) higher-order effects. It can be automatically computed by solving a modified 
version of Eq. (H), which is easier than performing perturbative calculations at 
high orders. Calculations of the neutron-deuteron (nd) phase shifts are presently 
available up to N^LO ( [T^ , |3^ ,|37D. At this order the only inputs are 7 and r^t (for 
the three-body quartet channels), and l/a^ and vqs (needed only for the doublet 
channels). There are no unknown LECs appearing at N^LO, so this approach 
can be easily pushed to higher orders (and precision). For a flavor of the results, 
we show the predicted quartet S'-wave phase shift in Fig. (|5|), compared to a PSA 
( p3| ) (at finite k) and a scattering length measurement (^^(essentially at /c = 0). 

Most of the data is not precise enough to provide a strict test for the con- 
vergence of the low-energy expansion, but the zero-energy point is much better 
measured. The EFT calculation gives for the quartet S'-wave scattering length 
03/2 = 5.09 + 0.89 + 0.35 -h. . . = 6.33 + 0.05 fm (P,|4|,|l2|), while the experimental 
value is 03^2 — 6-35 + 0.02 fm (|4^). The EFT error is probably overestimated 

since the N^LO calculation resummed some of the N'^LO pieces and the remain- 
ing ones are known to be smaller than the naive estimate (like the effect of the 
two-body shape parameter). Since the whole input in these calculations were 
the threshold parameters of A^A^ scattering, these results are universal (a "low- 
energy theorem"). Any model with the correct scattering lengths and effective 
ranges (and not wildly wrong phase shifts above threshold) should reproduce 
them, within the estimated error. The small discrepancy with "first-generation" 
NN potentials can be explained by the imprecise values of Cs^t that those models 
predicted. "Second-generation" (or "realistic") potentials fixed this discrepancy 

D- 

Coulomb interactions are important in proton-deuteron scattering at small 




Figure 6: Zero-energy half-ofF-shell amplitude for boson/two-boson scattering as 
a function of the outgoing momentum p, from a numerical solution of Eq. ( ]3l| ) 
with no three-body force, for several different cutoffs. 



energies. In Ref. (47) it was shown how they can be easily incorporated in the 



quartet channels by a simple change in the kernel of Eq. (^2|) . 
2.2.2 Triton-'^ He channel and the limit cycle 

The change of sign of A between the bosonic and quartet equations qualitatively 
changes the behavior of respective solutions. The most striking feature is that 
the solution of the bosonic (and of the S-wave doublet) equation in the absence of 
a three-body force depends sensitively on the value of cutoff used. This happens 
despite the fact that, as mentioned above, there is no ultraviolet divergence in 
the graphs summed. To illustrate this we show in Fig. (|^) the k = solutions of 
Eq. (|3l|) corresponding to various cutoffs (with X = 1, H = 0). Notice that the 
three-body scattering length, determined by the asymptotic plateau on the left, 
tf^=Q(k = 0), can take any value by just varying the cutoff within a reasonable 
range. 

It has been known for a long time that Eq. (^) is not well defined in the limit 
A — > oo. This disease can be shown in a variety of ways (|8|,|9|,|5|,|2l|). Consider 



the intermediate-momentum regime Q <C p <C A (remember that Q stands for 
any infrared scale Q j ^ k). Up to terms suppressed by powers of Q/A, and 
assuming ~ 1, Eq. (31) reduces to ( p8| ) 



°° 1 ( +pq + 



tk{p) = —f=- \ dq-\u \ tk{q). (33) 

V37rJo q \p^-pq + q^l 
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Eq. (|33[) has two symmetriesP 



ik{p) ti^{ap) (scale invariance) 

tk{p) (inversion symmetry), (34) 

that suggest the power-law solution tk{p) ~ p*. The allowed values of s are 
determined by plugging this ansatz back into Eq. ( |33[ ) . For values of s such that 
Re(s) < we find 

^ ^ - 1 (35) 

For values of A smaller than Ac = 3\/3/4'7r = 0.4135 . . ., the roots of Eq. ( |35| ) 
have Re(s) < —1 and the half-off-shell amplitude goes rapidly to zero as p — > oo. 
In the quartet channel, for instance, tk{p) ~ l/p^'^"^ , which is a softer ultraviolet 
behavior than the one expected in perturbation theory, t^ip) ~ For the 

S'-wave doublet case (or the bosonic case) though, A = 1 and there is a pair of 
imaginary solutions s = ibiso, with sq = 1.006. . . The asymptotic behavior of 
the half-off-shell amplitude is then 

tk{Q -^p <. A) = A sin{solnp + 6) (36) 

where A and 5 are p-independent numbers. This oscillatory behavior can be seen 
in Fig. d). 



The normalization A is fixed by the inhomogeneous term in Eq. (31). Because 
that is important only in the infrared region p ~ Q, j4 can be determined only 
by matching Eq. (|3^ ) to the solution in the infrared region. The phase 5 is 
determined by matching Eq. (^) to the solution in the ultraviolet region and 
will depend on the ultraviolet physics (like the value of the regulator A and the 
three-body force H{A)). Actually, by dimensional analysis, the A dependence of 
5 has to be the form (5(A) = —sq In A + 5. If A is varied while keeping H constant, 
as we have done so far, the dependence on the asymptotic phase 6 will "spill 
over" the infrared region and change the on-shell amplitude by a factor of 0(1). 
On the other hand, 5 does not depend on the infrared scales {k, 7, l/a^) and 
so H{A) can be adjusted in such a way that 6 = —sqIuA + 5{H{A)) = — sqIuA 
is cutoff independent (for any value of k), with A being a constant. Thus, the 
requirement of cutoff invariance means that H{A) runs at LO and, for generic 
values of A, H{A) ~ 1. The typical size of the three-body force is then 

(4vr)^ff (4^)^ 

as opposed to the naive estimate in Eq. (^). This means that Do is enhanced 
in the presence of fine-tuning in the two-body sector by a factor of (A/7)^. The 
arbitrary parameter A has to be fixed by one piece of experimental input coming 
from a three-body observable. The two-body phase shifts are not enough to fix 
the three-body physics already at LO. Another way of looking at the renormal- 
ization of the three-body system that is more easily generalizable to higher-order 
calculations is to consider two amplitudes t^{p) and (p) obtained by solving 

*These symmetries suggest that in the hmit 7, l/os — > the full conformal symmetry holds. 
In the two-body sector it was shown that this is indeed the case, even for off-shell amplitudes 
(|5l[). It is not known whether the three-body amplitude is conformal in this limit. The three- 
body forces break scale invariance at the order they first appear. 
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Eq. ( |3lD with two cutoffs A and A' and the corresponding three-body forces H{A) 
and H(A'). The integral equation determining (p) can be written as 

£ dq {5{p -q)- ^q^A{q)Kip, q)^ (g) 

2H{A) 22H{A) /"^ , o . ^ ^ A' ^ ^ 
= K{p, k) + — ^ + dqq^A{q)tt {q) 

- I dqq^^iq)tt (q) 

+ 1 l^'dqq'Aiq) g) + ^) (38) 

The first two hues in Eq. (^) are identical to the equation determining t^{p). 
If the effect of the remaining terms are small (up to terms of 0{Q/A)), t^{p) = 
t^' {p) (again, up to terms oi 0{Q/A)). These terms are indeed small, suppressed 
by a factor of Q/A compared to the leading ones. However, their effect on t^{p) 
is not suppressed. That is because the operator acting on t^{p) on the left-hand 
side of Eq. (|38| ) is nearly singular, that is, it has a small eigenvalue of order Q/A. 

£ (<5(p -q)- ^q^Aiq)Kip, q)^ a{q) = 0{Q/A)a{p) (39) 

The eigenfunction a{p) has the same asymptotics shown in Eq. (|36|) as t^{p). In 
our determination of the asymptotics, Eq. (|3^), we have already shown Eq. (|39| ) 
in the limit A — > oo and the total energy — > 0. We now see that the operator 
in the left-hand side of Eq. ( |38|) is almost non-invertible and that the projection 
of t^{p) in the a{p) direction is sensitive to the 0{Q/A) terms in the right side of 
Eq. (|38|). The amplitudes computed with different A's shown in Fig. (^) indeed 
differ in the intermediate regime by a term of the form sin(soln/> -|- 6). The 
solution t^{p) has an increased sensitivity to the ultraviolet physics and changes 
by a factor of 0{1) if the right-hand side of Eq. (^) changes by a factor of 
0{Q/ A). That explains the apparent contradiction between the cutoff sensitivity 
and the absence of ultraviolet divergences. 

Using the asymptotic form of t^{p), Eq. (p6|), even in the region p ~ A where 
it is unjustified, and dropping the terms suppressed by (Q/A)^ we can derive an 
approximate analytical expression for H{A) needed to cancel the Q/A terms in 
Eq. ( |38[ ) (and to guarantee t^{p) is cutoff independent up to order Q/A): 

cos (So In 4- + arctan sq ) 

H{A) = — ) \ (, (40) 

cos ( So In -j^ — arctan so ) 

where A is the parameter determining the asymptotic phase, to be fixed by ex- 
periment. We show H{A) in Fig. The points there were determined by 
numerically finding values of H{A) that, when inserted in Eq. (^), reproduce 
the same three-body phase shifts for different values of A. The solid line is 
Eq. (|40|) . The agreement between Eq. (|40| ) and the numerical points, as well as 
the independence of the phase shifts with A after H{A) is properly changed, con- 
firms our understanding of this somewhat unusual renormalization. For further 
discussion of renormalization-group invariance at this order, see Ref (|52[). 
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Figure 7: Three-body force coefficient H{A) computed analytically (blue line) 
and numerically (black dots) as a function of log (A/400 MeV). 

The asymptotic running of the three-body coupling can be interpreted as a limit 
cycle. The possibility of limit cycles in addition to fixed points was suggested in 
Ref. (^), but apparently never before seen in a simple physical system. Limit 
cycles are now being further studied (U). 

NLO calculations of phase shifts were performed in Ref. (|55| ) (scattering lengths 
were computed by different, but equivalent methods in Refs. (^ , |45|) ). No new 
three-body force is needed at this order (although H{A) has a correction pro- 
portional to the two-body effective range). As in other channels, it is easier to 
compute higher-order corrections by computing the kernel at the order desired 
and iterating it. That is, one numerically solves Eq. (^) with a kernel that 
includes higher-order effects. At NLO we demand t^{p) to be independent of 
A up to terms of 0{Q/A)'^. Due to the almost singular nature of the integral 
equation, the right-hand side of Eq. ( ^8|) has to be A independent up to terms of 
order {Q/Kf. This can be accomplished with the same no-derivative three-body 
force because the terms in the right-hand side of Eq. ( [3^ ) of 0{Q/K)^ are k, 
p independent. (However, H{A) will have a different form than in Eq. (^0[)). 
Similarly, at N^LO we need to cancel the 0{Q/A)^ terms in the right-hand side. 
Those terms, however, are proportional to k'^, and can be absorbed only by 
a three-body force with two derivatives. So, at N-^LO, a new three-body force 
appears, and we need yet another piece of three-body data to fix this new LEG. 
The cutoff sensitivity is dramatically reduced by going to higher orders, as it 
should (^6|). In Fig. (|8|) we show the phase shifts in the 5- wave doublet channel 
computed at the first three orders (|5^). As inputs in these calculations we have, 
besides the two-body interactions, a no-derivative three-body force fitted to the 
experimental binding energy of the triton at LO and a two-derivative three-body 
force fitted to the experimental value of the doublet nd scattering length at NLO. 
The EFT results are compared to a PSA (^) and to results from the Argonne 
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Figure 8: nd phase shift (in degrees) in the doublet 5-wave channel as function 
of the center-of-mass momentum: EFT at LO (light-blue dotted line), at NLO 
(blue dashed line), and at N^LO (purple solid line), a PSA (red crosses), and a 
potential model (black dots). From Ref. (|56|). 

V18 potential plus the Urbana three-body force adjusted to reproduce the correct 
triton binding energy (|57|). 

The existence, already at LO, of a parameter not determined by A^A^ scattering 
means that models tuned to reproduce the low-energy A^A phase shifts may differ 
widely on their predictions for three-body properties. However, since up to NLO 
there is only a one-parameter arbitrariness in the three-body predictions, there 
must be correlations among these predictions. This was noted empirically in Ref. 
(|58|). Fig. dH) shows results for the doublet 5- wave scattering length and the 
triton binding energy from a number of models, all reproducing the same low- 
energy A^A scattering (|59|) . The predictions cover a wide range but show a clear 
correlation ("Phillips line"). Also shown are the LO and NLO predictions for this 
correlation, obtained by varying the value of the three-body force at fixed cutoff. 
An equivalent explanation for the existence of the Phillips line was first found 
in Ref. (|60|). Analogous results can be obtained in the EFT for the hypertriton 

(H). 

In the 3A system, the pionless EFT seems to converge over a range of momenta 
that is large enough to include the interesting physics associated with the bound 
states. For example, if the three-body force is fitted to the scattering length, 
binding energy of triton is B3 = 8.08 + 0.23 + ... = (8.31 + . . .) MeV (||,|6|), to 
be compared to the experimental result B^^^ = 8.48 MeV. The success of these 
EFT calculations opens the way for the study of low-energy reactions involving 
triton and ^He. 

3 EFT WITH EXPLICIT PIONS 

As the typical momentum Q approaches the pion mass 771,^ it becomes increas- 
ingly difficult to account for pion exchange as a short-range effect. As we further 
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Figure 9: Correlation between the doublet 5-wave Nd scattering length and 
the triton binding energy (Phillips line): predictions of different models (black 
points), EFT in LO (light dashed line) and NLO (dark solid line), and experi- 
mental value (red cross). 



increase momenta past Q ~ mjr, we have to include in the EFT an explicit pion 
field and build up all its interactions allowed by symmetries. Because numerically 
the mass difference between the delta isobar and the nucleon, 6m = — rrifyf, 
is ~ 2m7r , convergence of the "pionful" EFT is optimized by the concomitant in- 
clusion of an explicit delta degree of freedom. The delta can be included without 
additional problems, since at these momenta the delta is, like the nucleon, a non- 
relativistic object. All other degrees of freedom can be considered heavy. Their 
effects are still subsumed in contact interactions, as they were in the pionless 
EFT. What we are doing is to remove part of the pion (and possibly the isobar) 
contributions from the contact interactions. One hopes the new EFT works for 
momenta up to Mqqd ~ 1 GeV, the mass scale of the heavier particles. 

Adding explicit pions to the theory will generate all sorts of non-analytic con- 
tributions to nuclear amplitudes. We want to devise a rationale for a controlled 
expansion in the presence of pions. 

3.1 Chiral symmetry and Chiral Perturbation Theory 

Fortunately, pion interactions are not arbitrary. Once explicit pion fields are 
considered, approximate chiral symmetry imposes important restrictions on the 
way pions couple among themselves and to other degrees of freedom. 

In the limit where we neglect the masses and m^ of the up and down 
quarks ("chiral limit"), QCD has a chiral SU{2)l x SU{2)r symmetry. Since 
this symmetry is not apparent in the QCD spectrum, it is reasonable to assume 
that it is broken spontaneously down to its diagonal subgroup, the SU {2)lj^r 
of isospin. Goldstone's theorem ( |6^ ) tells us that massless Goldstone bosons, 
naturally identified as pions, are associated with the three broken generators. 
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and their fields tt live on the "chiral circle" (actually a three-sphere S^) that 
represents the set of possible vacua. We call the radius of this sphere fj^; it is 
directly related to the chiral-symmetry-breaking scale, A^^SB ~ ^"^f-K ~ Mqcd, 
but the precise factor can only be obtained from the (so far elusive) mechanism 
of dynamical symmetry breaking in QCD. It turns out that this diameter can 
be determined from pion decay, and is called the pion decay constant, fj^ ~ 92.6 
MeV. 

The "chiral circle" can be parameterized in different ways but it is convenient 
to choose fields for which chiral symmetry is respected term-by-term in the effec- 
tive Lagrangian. Since the interactions of pions have to be invariant under chiral 
rotations, it is possible to choose fields where an infinitesimal rotation is repre- 
sented as TT ^ TT + e. Interactions of these fields always involve derivatives on 
the sphere, that is, derivatives together with certain non-linear self-interactions. 

As long as quark masses are small enough, their only effect is to change this 
picture slightly. A common quark mass breaks chiral symmetry explicitly down 
to the diagonal subgroup. Points on the chiral circle are no longer degenerate in 
energy, and a particular minimum is selected, in a direction given by the quark- 
mass terms that we define as the fourth direction. The quark mass difference 
further breaks isospin explicitly. In the low-energy EFT, the effect of quark-mass 
terms can be reproduced if we construct all terms that break chiral symmetry 
in the same way. These interactions can involve tt without derivatives, but are 
always accompanied by powers of rriu + rrid or niu — rrid. One example is a pion 
mass term, oc (m^ + m^). Likewise, electromagnetic and weak interactions 
can be constructed as well. 

The well-studied theory of non-linear representations of symmetries ( |63[ ) pro- 
vides the tools to write down the appropriate interactions between pions and 
other fields. The resulting chiral Lagrangian C has an infinite number of terms 
that can be grouped according to the index A: 

oo 

/:=^£(^), A = d + //2-2, (41) 

A=0 

where d is the number of derivatives, powers of m^r and/or powers of 5m, and 
/ is the number of fermion fields. Because of chiral symmetry, pion interactions 
stemming from QCD bring derivatives and/or powers of the pion mass. As a 
consequence, A > 0. The explicit form of for the lower values of A can be 
found in the literature (|6^ , |l^ ). 

As in the pionless EFT, the only hope of any predictive power rests on finding 
some ordering of contributions to amplitudes. This can be done for processes 
where all the external three- momenta are Q ~ m,r • Powers of Q of any particular 
Feynman diagram can be counted as done for the superficial degree of divergence. 
Each space derivative in an interaction produces a three-momentum in a vertex 
and therefore counts as Q. K complication in the counting of energies stems from 
the presence of heavy particles such as the nucleon together with light particles 
such as the pion. In any loop, integration over the zeroth component of the four- 
momentum will involve two types of poles, according to the scales appearing in the 
propagators: (i) standard poles at ~ Q corresponding to external three-momenta 
and to the mass of the pion; and (ii) shallow poles at ~ Q^/2m7v corresponding 
to external nucleon energies. 

Processes that involve at most one heavy particle {A = 0, 1) are the simplest 
because the contour of integration can always be closed so as to avoid shallow 



26 



Bedaque & van Kolck 



poles. In this case all energies are ~ As a consequence, each time derivative 
counts as Q and four- momentum integration brings a factor Q^. A pion prop- 
agator is and a nucleon (or delta) propagator is from its static term; 
kinetic terms are of relative 0{Q /rriN) and thus can be treated as corrections. 
With these ingredients one can write the contribution of any diagram to the 
amplitude as 

TocQ^J^iQ/A), (42) 

where A is a renormalization scale, .7-" is a calculable function of LECs, and ly is 
a counting index. For strong interactions (|65|), 

iy = A-2C-A + 2L + ^Ai, (43) 

i 

where C = 1 is the number of connected pieces, L is the number of loops, and the 
sum runs over all vertices. In addition, electroweak interactions can be considered 
through a simultaneous expansion in a = e^/Air and Gpf^. Since L > and 
A>0, i'>2 — A = Vmin- Assuming that all LECs have "natural" size (given 
by powers of Mqcd once the lower scales have been identified explicitly), an 
expansion in Q/Mqcd results. Its first two orders are equivalent to the current 
algebra of the 60's, but at higher orders unitarity corrections can be accounted for 
systematically. In the sector of j4 = 0, 1, the EFT is called Chiral Perturbation 
Theory (ChPT). 

In processes that involve more than one stable heavy particle {A > 2), on the 
other hand, a failure of perturbation theory can lead to bound states (|6|). The 
shallow poles cannot be avoided: they represent "reducible" intermediate states 
that, as in the pionless EFT, differ from initial states only by nucleon kinetic 
energies, which are ~ Q^/m-Tv- This 0{mj\[/Q) infrared enhancement over inter- 
mediate states where energies are ~ Q invalidates Eq. (E^) for reducible diagrams. 



Fermion lines, loops and derivatives then scale with Q as in Eq. (16). A pion 
propagator still counts as Q~'^, but the pion can be taken in first approximation 
as static, and it is sometimes referred to as a "potential" pion. Contributions 
that come from standard poles naively scale as in processes with no more than 
one heavy particle. Pions there are non-static or "radiative" . 

The issue now is how to estimate the size of pion contributions. One needs 
to find the importance oi (i) pion exchange relative to short-range interactions; 
(ii) multi-pion exchange relative to one-pion exchange (OPE). Both issues are 
related, via renormalization, to the size of the contact interactions. How large 
are they in the pionful EFT? What are the contributions that must be resummed 
in nuclear amplitudes? 



3.2 The two-nucleon system 

As we have discussed, the A^A^ system is characterized by scattering lengths ag, 
at that are much larger than 1/Mqcd- In the pionless theory, this fine-tuning 
cannot be explained, but it can be accommodated in the power counting by 
assigning to the contact interactions the scaling given in Eq. (|l^). 

A new scale appears naturally in the pionful theory. The leading (A = 0) 
coupling of the pion to the nucleon is derivative with a coupling constant qa/ fn, 
where — 1.26 is a parameter not fixed by symmetry but determined in /3-decay. 
The OPE contribution to the NN amplitude is, schematically, g'j^Q'^ / /^[Q'^+m'^). 
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Because a reducible intermediate state contributes mjyQ/'iTr, relative to OPE, 
once-iterated OPE can be estimated to give a contribution 

9\QymQ^ + ml) ^ Mnn ^ ' 

(as long as Q ^ w-tt)- Here we introduced the scale 

Mnn = (45) 

which sets the relative strength of multi-pion exchange. Numerically (for A^^^ = 3) 
~ 47r/7r and ~ 1, so M^n ~ /tt- This naive dimensional analysis cannot, 
however, capture the numerical factors that actually set the relative size of pion 
contributions. A more accurate estimate requires concrete calculations. 



3.2.1 Perturbative pions 

For Q < M]\fN, iteration of OPE should be suppressed with respect to OPE 
according to Eq. (^^. Moreover, if we assume that the leading short-range 
effects are ~ Aira/rriN (as in Eq. ( p^ with Q ^ 1/a), then OPE is suppressed by 
0{l/aM]\fj\f) compared to the leading contact interaction. With such estimates, 
if Mat TV is sufficiently large (compared to 1/a and 771,7^) and Q sufficiently small, 
pions may be treated perturbatively. The suggestion that this can be profitably 
done was made in Refs. (|66|, p^). 

A simple power counting, which became known as KSW counting, follows by 
taking Q ~ 1/a ~ m-,^ < Mjy^, and counting powers of the light scale Q. This is 
a direct extension of the power counting in Eq. (0). In particular, the scaling 
of the contact operators is assumed to be the same as before with A M^n, 
and thus their ordering is unchanged. Because of chiral symmetry, each insertion 
of a pion exchange brings a factor of Q/M^n- Electroweak interactions can be 
treated in much the same way as before. One can show that renormalization can 
be carried out consistently within this power counting (p3|). 

However, for this power counting to be relevant to nuclear physics Mnn has to 
be sufficiently large. If Mtvat is not larger than ttItt, the domain of perturbative 
pions is no larger than that of the simpler pionless theory. The issue of the 
range of validity of the EFT with perturbative pions can only be settled by 
explicit calculation of dimensionless factors and comparison with precise observed 
quantities, this being done to sufficiently high order so that a significant number 
of pion effects be tested. 

With this power counting, the LO A^A^ amplitude coincides with that in the 
pionless EFT, see Fig. (|l|). At this order there are contributions only in the two S- 
wave channels from chir ally-symmetric, non-derivative contact interactions (the 
Cq terms). Subleading terms, of 0{Q/MN]y) relative to leading, are constructed 
in a direct extension of subleading terms of the pionless EFT. Besides two two- 
derivative contact interactions (C2 terms), we also insert OPE and two non- 
derivative contact interactions that break chiral symmetry explicitly {m^Cf^ 
terms). Both C2 and Cf^ contact interactions only contribute to S waves. The 
tensor operator from OPE, on the other hand, introduces mixing between ^5*1 
and waves. To this order, all but the S'-wave phases are predicted in terms 
of pion parameters. A calculation of the A^A^ system to NLO was carried out in 
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Figure 10: ^Si NN phase shift in the EFT with perturbative pions, as function 
of the center-of-mass momentum. The long-dashed, short-dashed, and dotted 
lines are, respectively, the EFT results in LO, NLO, and N^LO. The dash-dotted 
line is the N^LO result with a higher-order contact interaction added. The solid 
line is the Nijmegen PSA. From Ref. ([72|), courtesy of T. Mehen. 



Refs. (^,67). When the results were compared to the Nijmegen PSA (|2q), it 
was suggested that M^r n ~ 300 MeV (p3| . The relative size of OPE and contact 
interactions was extensively discussed (|68,^,70), but the issue is clouded by the 
details of fitting procedures. 

At N^LO we encounter new pion exchanges: both non-static (or radiative) 
OPE and once-iterated OPE. A calculation at this order was carried out in all S, 
P, and D waves (72,71). Results were found to depend on the channel. While in 
singlet channels, like ^Sq, there seems to be good agreement with the Nijmegen 
PSA, in triplet channels, such as ^Si, ^Di, and ^^0,2, the N'^LO corrections are 
big at momenta ~ 100 MeV and lead to large disagreement. One example is 
shown in Fig. (p!o|). The effects of perturbative pions (and of delta isobars) is 
milder in the higher partial waves (73). 

The problem can be traced to the iteration of the tensor part of OPE. These 
results suggest that pions, or more explicitly, the Yukawa part of potential and 
radiation pions, when treated perturbatively give rise to a converging expansion 
for the ^5*0 scattering amplitude up to fairly large momenta. However, OPE in 
the ^Si — ^Di coupled channels is not perturbatively convergent for momenta 
around 100 MeV, because the tensor force, which survives in the chiral limit, 
is too large. This, in turn, suggests that the naive estimate Mnn ~ fn is not 
entirely unreasonable. 
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3.2.2 Renormalization of the pion ladder and power counting 

If indeed Mattv <Si Mqcd, we might hope to improve on the expansion of the 
previous section by a controlled resummation of terms that go as Q/M]\jn. If 
Mnn ~ 100 MeV, the lack of other known particle thresholds there suggests that 
the resummation could involve only pions. Indeed, Weinberg, who first attempted 
the use of EFT in the derivation of nuclear forces (|^, had already suggested a 
power counting in which pions appear in leading order, and should therefore be 
iterated. (Some elements of this power counting had been anticipated by Priar 
(0))- 

Weinberg's original proposal (^ for an EFT describing multi-nucleon systems 
was to split the full amplitudes into reducible and irreducible parts. Irreducible 
diagrams, in which typical energies superficially resemble those in ordinary ChPT, 



should satisfy the power counting, Eq. (43). We call the sum of irreducible 
diagrams the potential V . Note that the potential, being a set of subgraphs, can 
be defined in alternative ways. All definitions that exclude the infrared-enhanced 
contributions but differ by a smaller amount are equally good, as long as no 
double-counting or omissions are made. A field redefinition might change the 
potential but not the full amplitude. The important point is that the only scale 
appearing explicitly in the potential is Q, so that the power counting proceeds as 
in the case of diagrams with at most one heavy particle. Reducible diagrams can 
be obtained by sewing together irreducible diagrams via intermediate states that 
contains the propagation of only the initial particles. The full amplitude T for an 
^-nucleon system is thus a sum of the potential and its iterations; schematically, 

T = V + VGoV + VGqVGqV + ... = V + VGqT, (46) 

where Go is the ^-nucleon free (Schrodinger) Green's function. This is just the 
Lippmann-Schwinger equation, which is formally equivalent to the Schrodinger 
equation with the potential V, from which wave-functions can be constructed. 

Unfortunately, just counting powers of Q is not in itself sufficient for an ordering 
of interactions. We need to find which other scales accompany Q. Given that we 
do not yet possess a full solution of QCD, it is obvious that some assumptions 
have to be made about the LECs. Because OPE has short-range components, it 
is natural to assume that they set the scale for the short-range interactions. In 
this case, since OPE is 0{g\/ f"^) for Q ~ m-,^, we expect the leading two-nucleon 
contact interactions to be 0{4:TT/m]\fMN]\f). This assumption naturally explains 
why nuclear bound states are much shallower than naively expected (HP). The 



series in Eq. (46) is roughly 

„ 47r 



uinMnn 



(47) 



which requires resummation and exhibits a (real or virtual) bound state at Q ~ 
MjvTv- In other words, the natural scales for the A^A^ scattering length and for 
the binding energy of a nucleus are 

H~T7^, B^^, (48) 
Mnn ruN 

respectively. If we use again that miy ~ 'i'^fn and qa ~ 1, then B ~ f-w/^T^ ^ 10 
MeV. So we find that it is not Mqcd by itself that sets the scale for binding 
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energies, but a certain ratio of powers of f-j^ and niN. Now, as we have seen, this 
is not the whole story. It remains mysterious why the A^A^ (real and virtual) 
bound states are even shallower, or equivalently, why \a\ > X/M^n, by factors of 
a few. This still has to be accommodated by fine-tuning the contact interactions. 

Because they serve as counterterms to pion loops in the potential, which are 
expected to be suppressed (as in ChPT) by powers of Q/^T^f-K, Weinberg implic- 
itly assumed that LECs related to more derivatives and powers of the pion mass 
contain inverse powers of Mqqd- That is, a (renormalized) contact operator with 
index A would scale as 

Ca . , (49) 

as in naive dimensional analysis. 

With this assumption and disregarding the fine-tuning, a simple power count- 
ing results from taking Q ~ M]\fj\f ~ 1/a ~ m^. The potential obeys Eqs. 
(^,43). The leading potential consists of no-derivative, chirally-symmetric con- 
tact interactions plus static OPE, 



where ai (2tj) are the Pauli matrices in spin (isospin) space and qis the momen- 
tum transferred. All contributions to nuclear forces other than Eq. ( |50[ ) would 
come as corrections in powers of Q/Mqcd, starting at {Q /MqcdY ■ The struc- 
ture of the potential rapidly becomes more complex with increasing order (^. 
The leading potential has to be resummed, while corrections can be treated in 
perturbation theory. If the corrections are truly small, resumming them should 
cause no major harm. This method requires numerical solution of the Schrodinger 
equation and is similar in spirit to the traditional potential-model approach. As 
we will see below, Weinberg's power counting has been extensively and success- 
fully developed during the past decade to study processes involving few-nucleon 
systems. 

However, there is a subtlety not present in ChPT. Loops in reducible diagrams 
probe, as other loops do, high energies (when nucleons are far off shell). As in the 
pionless EFT, the potential does not vanish at large momenta: it is singular. In 
addition to the delta function and its derivatives already present in the pionless 
theory, pion exchange generates potentials that behave as l/r" with n > 3, as the 
radial coordinate r — > 0. The large-momentum or short-distance behavior is the 
same as in the chiral limit 0. Already in leading order the tensor force goes 

as in the chiral limit. As a consequence of ultraviolet divergences generated 
by the iteration of the potential, the infrared enhancement of Mqcd/^nn might 
contaminate LECs, possibly invalidating Eq. (^). 

The crucial issue is whether, at any given order, all divergences generated by 
iteration can be absorbed in the parameters of the potential truncated at that 
order. There is some indication that equally good fits can be achieved in leading 
orders with various cutoffs (^,|75,76, 70|,^,|7^), as required of a sensible EFT, but 



the numerical nature of the results makes a definite answer difficult. 

Unfortunately, there seem to be formal inconsistencies in Weinberg's counting 
( p3| , |69| ) . Divergences that arise in the iteration of leading-order interactions 
apparently cannot be absorbed by the leading-order operators themselves. Two 
examples are two-loop diagrams where: i) OPE happens between two contact 
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interactions, having a divergence proportional to the square of the pion mass (|2^ 
ii) OPE is iterated three times, having a divergence proportional to the square 
of the momentum (^9|). Although these two particular cases could be resolved 
by the promotion to leading order of two counterterms, it is likely that a similar 
problem would show up at higher orders in the expansion. The correspondence 
between divergences and counterterms appears to be lost, a fundamental problem 
with the chiral expansion and the momentum expansion resulting from Weinberg 
power counting. 

This argument is not decisive, though. It has been known that in the context 
of the Schrodinger equation perturbative arguments are not in general reliable 
for singular potentials (^). The perturbative expansion might have a cut start- 
ing at g\/ f"^ = 0; insistence on a g\/ f"^ expansion would then reflect itself on 
different orders offering correlated contributions to counterterms, each bringing 
powers of Mqcd/M^nn yet resulting in a much better behaved sum. How can the 
resummation be consistent? There is a mapping between the singular two-body 
1/r^ central potential and the three-body problem with short-range interactions. 
In Sect. ^J2?\ we saw that for the latter, the renormalization of the nonpertur- 
bative equation is very different from the renormalization of individual terms in 
the associated perturbative series (|2l|). In particular, the relevant counterterm 
exhibits a limit-cycle behavior. 

It turns out that the correct renormalization of singular potentials is, indeed. 



in general intrinsically nonperturbative (^, 8C). In contrast to regular potentials, 
both solutions of the Schrodinger equation for an atractive 1/r" central potential 
are equally acceptable: the radial wave- function n(r), whose large-distance be- 
havior determines low-energy observables, oscillates rapidly as r decreases. Per- 
turbative approximations to the wave-function fail at distances comparable to 
the intrinsic scale ro present in the potential, as illustrated in Fig. (11). The 
problem can be rendered essentially cutoff independent with a single counterterm 
associated with the short-distance physics. For example, short-distance physics 
can be represented by a square well potential of radius i? <C whose depth 
Vb = yo{R) can be adjusted so that physics at r ^ tq be independent of R (|80|) . 
(The advantage of this coordinate-space regulator is that one can do an analytic 
matching of the outer and inner solutions of the Schrodinger equation, thus find- 
ing the desired Vq{R). For another technique to deal with a singular potential, 
see, for example, Ref. (|8l|).) While for a repulsive potential there exist only 
fixed points (^2|), the situation in the attractive case is similar to the three-body 
system. 

These results were extended to the A^A^ potential in the ^Sq and — ^Di 
channels in Ref. (|8^), where the long-range part of the potential was taken as 
that of OPE for r > R. The asymptotic behavior is that of the chiral limit, 
where the relevant scale is tq ~ 1/Mnn- The depths of the short-range part of 
the potential can be different in the singlet and triplet channels, as there are two 
parameters (Cq*^^ and Cq"^^) in Eq. (]50|). 

In the ^5*0 channel the calculation is straightforward. The pion potential is 
simply a Yukawa, and the explicit solution is 

V,AR, n) = -(2„ + 1)^^ - ^ (I) + 0(B!'), (51) 

where i?* is an intrinsic length scale to be determined numerically from a fit to 
low-energy data, and n labels the branch of a cotangent. In the left panel of Fig. 
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Figure 11: Zero-energy radial wave-functions of the 1/r^ potential: exact (solid 
line) and perturbation theory to LO (small dashes), NLO (medium dashes), and 
N^LO (large dashes). Prom Ref. (^. 



(|T^ the R dependence of Vq, as given by Eq. (|5T|), is compared to the numerical 
solution of the Schrodinger equation with the observed singlet scattering length. 
The presence of a multi-branch structure is related to the accumulation of bound 
states inside the square well. Of course the presence of unphysical bound states is 
innocuous as long as the binding energies of such states are near the cutoff of the 
EFT. One sees that the formal problem with the chiral expansion in Weinberg's 
counting survives the resummation. While the first cutoff-dependent term in Eq. 
(51) can be represented by a chiral-symmetric contact interaction, the second 
would require a chiral-breaking one. In momentum-space notation where the 
cutoff is denoted by A, 



Co(A) 



Air 



6mjvA 



vr 



(2n + l)V + 



1 



(52) 



Although the logarithmic divergence is suppressed by a power of A compared 
to the first term, it is a true divergence in physical quantities that must be 
renormalized at leading order in Weinberg power counting. The C^^^ operator, 
which is formally subleading, must be promoted to leading order if the full OPE 
is iterated, in agreement with the perturbative argument of Ref. ([2^). On the 
other hand, the C^^'^ contribution is numerically small, as demonstrated by the 
dotted curve in the left panel of Fig. (jl^) which neglects the A~^ corrections to the 
running. This smallness explains why Weinberg's power counting has been found 
to work well in this channel over a moderate range of cutoffs (1,^^^^^. 

A possible conclusion is that OPE and the O {ruq) LECs contribute to any 
amplitude at the same order in the expansion, and this is what leads to KSW 
power counting. However, a more general conclusion to draw is that the difference 
between the OPE contribution for rUq ^ and the OPE contribution in the chiral 
limit must occur at the same order as the O {mq) counterterms. In many cases 
these two conclusions yield identical amplitudes, however, in the ^Si — ^Di channel 
they do not. 

In the ^S*! — channel, in addition to the long-distance Yukawa interaction 
and the contact interaction, there is a strong tensor component of OPE that cou- 
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R (fm) R (fm) 

Figure 12: Running of mj\fVQ as a function of the cutoff R. Singlet channel 
(left panel): Eq. (^) for n = 1 (blue solid line), and same but with part 
neglected (red dotted line). Triplet channel (right panel): numerical solution of 
an analytic matching equation (solid lines). Dots are extracted directly from a 
numerical solution of the Schrodinger equation in the respective channel. From 
Ref. (g). 



pies S and D waves. At distances r <C l/m,r the central potential is negligible, 
while in the region R < r <^ I/Mtvat we can neglect the angular- momentum 
barrier. Moreover, for \pmf^ <C Mat at the total energy E can be treated as a 
perturbation. In this short-distance limit, we can keep only the chiral limit of 
the tensor force, and the Schrodinger equation can be diagonalized and solved 
exactly. In the diagonal basis the Schrodinger equation decouples into an attrac- 
tive singular potential and a repulsive potential. The solution for the attractive 



singular potential is a linear combination of Bessel functions (34, 30), and the 
wave-function at this order is 



„M = ^.-»'^cos(^-^ + ^„j, (53) 

where ^ is a dimensionful normalization constant and </>o is the asymptotic phase 
which determines the triplet scattering length. This solution oscillates ever faster 
as it approaches the origin, just as in Fig. (lTl|). As before, the issue is whether 
a Vb(i?) can be found in such a way that the asymptotic phase is made R 
independent. Matching logarithmic derivatives of the interior square- well and 
exterior attractive solutions at r = 72 yields an equation whose solution is shown 
in the right panel of Fig. ([l2|) . The renormalization-group flow is multibranched 
and non-analytic in g\/ f"^. 

The leading-order ^Si phase shift is clearly cutoff independent as shown in 
Fig. (|l|). The situation is similar for the ^Di phase shift, while the mixing 
parameter ei exhibits some R dependence. However, an error plot of si suggests 
that the R dependence and the deviations from the Nijmegen PSA (26) are 



higher order in the momentum expansion. These results (|83D are in agreement 



with the numerical analyses of Refs. (34, 77). In these channels, contrary to 



the perturbative argument (p9D, Weinberg's power counting does not seem to be 
formally inconsistent. 

It was conjectured in Ref. (^3|) that a formally consistent expansion in the 
pionful EFT is an expansion around the chiral limit. This expansion is equivalent 
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Figure 13: '^Si NN phase shift in leading order in the EFT with non- 
perturbative pions, as function of the center-of-mass momentum. The (pur- 
ple) long-dashed, (red) medium-dashed, and (red) short-dashed lines have, re- 
spectively, R = 0.45 fm (A = 438 MeV), R = 0.21 fm (A = 938 MeV), and 
i? = 0.10 fm (A = 1970 MeV). The (black) sohd line is the Nijmegen PSA. From 
Ref. (g). 



to KSW power counting in the channel and equivalent to Weinberg power 
counting in the ^Si — ^Di coupled channels, i.e. it selects only the desirable 
features of both power countings. The leading-order potential y(r;0), to be 
treated exactly, consists of the chirally-symmetric component of OPE and non- 
derivative contact interactions. Deviations from the chiral limit V{r] m^)— y(r; 0) 
can be treated perturbatively in all channels, and in fact, such a perturbative 
expansion is required in the ^5*0 channel but not in the ^Si—^Di channel. Evidence 
was presented that this expansion converges, even though convergence is slow 
due not to the long-range pion physics itself, but to the fine-tuned short-distance 
physics (as argued previously (|85|)). 

Although existing obstacles to a derivative and pion-mass expansion were re- 
moved, higher orders must be studied before the issue can be considered set- 
tled. For example, in Ref. (^) an incomplete sub- leading calculation with non- 
perturbative pions has found limits in fitting the effective range. Note also that 
alternative views of the renormalization of the pion ladder exist (p7|). Finally, 
there is an interesting suggestion of expanding the A^A^ amplitude in the en- 
ergy region where the S-wave phase shifts vanish (|88D. The connection between 
this expansion and the low-energy expansion described here has not been fully 
analyzed. 

3.2.3 Potentials and fits to A^A^ data 

The picture that emerges from the previous section is close, conceptually and 
numerically, to Weinberg's original proposal. To be formally consistent, we should 
expand in the pion mass. Yet, if we resum the effects of the pion mass in pion 
exchange, which can be done with higher-order error, the leading-order potential 
becomes the same as in Weinberg's power counting. Corrections to the leading 
potential do not need to be iterated to all orders. Yet, as has been shown explicitly 
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in the pionless EFT (0), they can be iterated with small error, as long as one 
uses a regularization with a cutoff A ~ Mqcd- Clearly, a potential which is 
correct up to a certain order ensures that the amplitude is correct to the same 
order. 

Much work has been done in developing an EFT potential based on Weinberg's 
power counting. Traditionally, potential models have been plagued by problems 
of principle, such as the form of meson-nucleon interactions (for example pseu- 
doscalar vs. pseudovector pion coupling), renormalization issues, absence of a 
small expansion parameter, etc. Because the EFT potential includes explicitly 
the exchange of only pions, all these problems can be resolved. For any given 
choice of pion field, the form of interactions is fixed by the pattern of chiral 
symmetry breaking. Renormalization can be performed because all interactions 
consistent with symmetries are included. And the power counting Eq. (43) for 
the EFT potential implies that diagrams with an increasing number of loops L 
— and, in particular, with increasing number of exchanged pions — should be pro- 
gressively less important. The EFT potential can be thought of as a low-energy 
approximation to standard potential models, although this can only be taken in 
an average sense. 

In leading order, u = Vrnin = 0, the A^A^ potential is, as we saw in the previous 
section, simply static OPE and momentum-independent contact terms. This is 
obviously a very crude approximation to the A"A^ potential: it is known that the 
nuclear force has other sizable components, like a spin-orbit force, a strong short- 
range repulsion and an intermediate range attraction. These are all generated in 
the next orders: v = Vmin + 1 corrections vanish due to parity and time-reversal 
invariance, but v = ymin + 2 corrections are several. First, there are short-range 
corrections; they come from one-loop pion dressing of the lowest-order contact 
interactions, and from four-nucleon contact interactions with two derivatives or 
two powers of the pion mass. It is easy to show that the result of loop diagrams 
amount to a simple shift of the contact parameters. Second, there are corrections 
to OPE; these come from vertex dressing and from recoil upon pion emission. 
Third, there are two-pion exchange (TPE) diagrams built out of lowest-order 
ttNN (and ttNA) interaction. At = I'min + 3 a few more TPE diagrams 
appear, which involve the ttttNN seagull vertices from the A = 1 Lagrangian. 
To this order there are also small some relativistic corrections. At = i'min + 4 a 
host of two-loop diagrams and new contact interactions emerge, and so on. Some 
diagrams are shown in Fig. (ll^). 

A calculation of all contributions to the A^A^ potential up to z/ = fmin + 3 was 
carried out in Refs. (0, |9|) using time-ordered perturbation theory. This EFT 
potential is energy dependent, but equivalent potentials can be obtained through 
unitary transformations. An energy-independent potential is more convenient in 
many situations, and the corresponding version was derived in Ref. (p9|). 

The potential to this order has all the spin-isospin structure of phenomenolog- 
ical models, but its profile is determined by explicit degrees of freedom, symme- 
tries, and power counting. The power counting suggests a hierarchy of short-range 

(S T) 

effects: S waves should depend strongly on the short-range parameters Cg ; 
contact interactions affect P-wave phase shifts only in subleading order, so their 
effect should be smaller and approximately linear; and D and higher waves are 
directly affected by contact interactions at higher orders, being thus essentially 
determined by pion exchange. While phenomenological potentials such as that in 
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Figure 14: Some time-ordered diagrams contributing to the NN potential in the 
pionful EFT. (Double) solid lines represent nucleons (and/or deltas), dashed lines 
pions, a heavy dot an interaction in C^'^\ a dot within a circle an interaction in 
C^^\ and a dot within two circles an interaction in C^'^\ The first line corresponds 
to 

— ^min^ second and third lines to v — I'min + 2, fourth line to — z^mm + 3, 
and ". . ." denote v > Vmin + 4. All orderings with at least one pion (or delta) in 
intermediate states are included. Not shown are diagrams contributing only to 
renormalization of parameters. 



Ref. ( |90D have similar short- and long-range structure, it is on TPE that chiral 
symmetry is particular influential. TPE here includes a particular form of terms 
previously considered (pll), plus a few new terms. These new terms come from 
nucleon-structure properties, such as the axial polarizability and the sigma term, 
and they provide a correlation between the exchanged pions that is important 
in the isoscalar central force. (Even though graphs where pions interact in flight 
appear only in next order and should thus be relatively small.) Not surprisingly, 
in the chiral limit these components of the potential behave at large separations 
as van der Waals forces. The components of the TPE potential were studied 
in detail in Ref. ([73|). In particular, it is shown explicitly that i) relativistic 
corrections are mostly small; ii) both isoscalar central and spin-orbit potentials 
are numerically similar to a and uj exchange in models; in) the OPE isovector 
tensor potential is reduced by the TPE contribution. More recently, (the tail 
of the energy-independent version of) the EFT TPE potential in the limit of a 
heavy delta was substituted for (the tail of) the one-boson exchange in a Ni- 
jmegen phase-shift reanalysis of pp data below 350 MeV (92). A drop in was 
observed. When rnrNN seagull LECs are fitted to the A^A^ data, they come out 
close to values extracted from ttN scattering. This confirms unequivocally the 
validity of chiral TPE. A new full Nijmegen PSA is in the works, in which chiral 
TPE is used in the long-range potential (|92|). 

Fits to A^A^ phase shifts were done to this order with (^) and without (|78| ) 
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Figure 15: "^Si NN phase shift (in degrees) in the EFT with Weinberg counting, 
as function of the lab energy (in GeV). The (purple) dotted, (green) dashed and 
(red) solid lines are the results at orders = i^min, ^ = + and u = fmm + 3, 
respectively. The squares are the Nijmegen PSA. From Ref. (j7^), courtesy of 
U.-G. MeiBner. 



explicit delta degrees of freedom. For numerical convenience, smooth cutoffs are 
used to regulate the loops generated by the Schrodinger equation. For each cutoff 
value a set of nine independent (bare) parameters stemming from contact inter- 
actions were fitted to phase shifts below 100 MeV laboratory energies. However, 
the fits differed in the details. In Ref. (|^, a Gaussian cutoff was used in all 
loops, and calculations performed with the cutoff parameter A taking values 500, 
780 and 1000 MeV. In contrast, Ref. (^) used a mixed regularization: while the 
cutoff parameter in the smooth regulator was varied between 500 and 600 MeV, 
loops in the potential were treated by a subtraction procedure equivalent to DR. 
Although OPE and TPE diagrams are completely determined by LECs accessible 
in TiN reactions, most of these LECs were not known at the time of Ref. (^), 
and were also searched in the fit. In Ref. (^), these parameters were taken 
from a (DR) fit to ttN scattering, which allowed for a simpler, partial- wave- by- 
partial-wave fit of A^A^ phase shifts. Reasonable agreement with existing PSAs 
and deuteron properties was found, especially in Ref. (|7q). As an example, the 



^Si phase shift at various orders is shown in Fig. (15) and compared to the Ni- 
jmegen PSA (^6|). (Cf. Figs. (PI, p!o|, p^.) The contributions from the short-range 
parameters in this fit turn out to be comparable to those from heavier resonances 
in phenomenological models (93). 

Although the fits to = Vmin + 3 are good, they are inferior to the so-called 
"realistic" potentials that use 40 — 50 parameters to fit data up to 300 MeV with a 
near 1. First steps are being made to extend the EFT potential to = fmm + 
4. In Ref. (p^), contact interactions of this order were added to the complete 
deltaless u = I'min + 3 potential, bringing the number of adjustable parameters 
to the level of realistic potentials and the Nijmegen PSA. The resulting fit is of 
quality comparable to realistic potentials. (Indeed, the difference among realistic 
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potentials arise from high-momentum intermediate states that probe scales which 
cannot be uniquely fixed by low-energy fits (|9^).) In order to achieve a complete 
1^ = z^mm + 4 potential, two-loop diagrams and corrections to one-loop diagrams 
need to be included. Classes of these diagrams that are invariant under pion-field 
redefinitions have been calculated in Ref. (96). They are compared to peripheral 



partial waves in Ref. (97), where they are found to be relatively small. We expect 



to soon witness a complete v = Vmin -|- 4 fit of high accuracy. 
3.2.4 Isospin violation 

The mass difference between u and d quarks breaks isospin symmetry explicitly. 
Indication from the meson masses is that the ratio e = (m^ — rrid) / {rriu + rrid) ~ 
1/3. Naively, this suggests that isospin might not be a much better symmetry 
than the rest of the chiral group. On the other hand, a cursory look at hadron 
masses and a more complete analysis of dynamical amplitudes show that isospin 
is typically broken only at the few percent level. 

Why is isospin such a good symmetry at low energies? The answer can be found 
in the pattern of chiral-symmetry breaking incorporated in the chiral Lagrangian 



(11). While explicit chiral-symmetry-breaking effects are present already at index 
A = through the pion mass term, operators generated by the quark mass 
difference appear only at A = 1 through a term that contributes to the nucleon 
mass splitting and, due to chiral symmetry, to certain pion-nucleon interactions. 
As a consequence, in most quantities isospin breaking competes with isospin- 
conserving operators of lower order, and its relative size is not e but £{Q/Mqcd)^, 
where n is a positive integer. In other words, isospin is an accidental symmetry 
([ll|) : a symmetry of the lowest order EFT which is not a symmetry of the 
underlying theory. The only known exception to this rule is in the isoscalar t 
channel in vrA'^ scattering at threshold. There, there is no contribution from the 
A = Lagrangian, and both the isospin-conserving and -violating amplitudes 
start at the same order. The isospin-violating piece comes from the pion-nucleon 
interactions linked to the nucleon mass splitting. Unfortunately, this is hard to 
see experimentally. 

Along these lines, one can study the expected size of isospin breaking in the 
nuclear potential. We follow the standard nomenclature and refer to an isospin- 
symmetric potential as "class I", to a potential that breaks charge dependence 
but not charge symmetry — defined as a rotation of vr around the 2-axis in isospin 
space — as "class II" , to one that breaks charge symmetry but vanishes in the np 
system as "class III" , and to one that breaks charge symmetry but causes mixing 
in the np system as "class IV". 

At (5 ~ MjsfN, photon exchanges are perturbative. These standard electromag- 
netic effects from "soft" photons can be obtained straightforwardly from gauge- 
invariant operators involving the photon field. In addition, isospin violation arises 
from the quark masses, from indirect electromagnetic effects, and from simulta- 
neous pion-photon exchange. In order to compare the various sources of isospin 
breaking, we note that the size of electromagnetic effects in loops is typically 
~ a/vr which, numerically, is ~ ^{Q /MqcdT ■ 

The leading isospin-breaking interactions in Weinberg's power counting have 
been derived in Ref. (11). (The necessary modification of Weinberg's counting 

shows that chiral-symmetry-breaking terms are even 



discussed in Sect. 3.2.2 



more suppressed vis-a-vis chiral-symmetric ones, and is not expected to affect the 
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relative sizes among isospin-breaking interactions.) No isospin- violating effects 
appear at leading order, u = Umim so the leading potential is class I. The first 
isospin-breaking effect (in addition to Coulomb exchange) appears at u = Vmin + 1 
in the form of a class II isospin violation from the pion mass splitting (Am^ = 
0{aMQ(jj-)/-K)) in OPE. One order down, v = Vmin + 2, a class III force appears 
mainly from the quark mass difference through breaking in the ttNN coupling 
{Pi = 0{eml/M^^^)) in OPE , from contact terms (7^,,^ = 0{eml/M^cD))^ 
and from the nucleon mass difference {Attin = 0{em'^/MQCD))- To this order 
the isospin-violating nuclear potential is a two-nucleon potential of the form 



where 



Vib = Vii [(^1)3(^2)3 - ti ■ <2] + ^iii [{ti)3 + (^2)3] , (54) 



Vn = -( ^V ,^,^%\7^;: 2 . (^rni + Am%), (55) 



= T^T^T^ " + ^-"^ ■ ^''^ 

Finally, class IV forces appear only at order v = Vmin + 3. 

We conclude that the pattern of symmetry breaking in QCD naturally suggests 
a hierarchy of classes in the nuclear potential ([ll|): 

(•^M+i) n(-9—\ 

where (^m) denotes the average contribution of the leading class M potential. 
This qualitatively explains not only why isospin is a good approximate symmetry 
at low energies, but also why charge symmetry is an even better symmetry. 
It gives, for example, the observed isospin structure of the Coulomb-corrected 
scattering lengths (||), a„p ~ 4 x {{cnn + app)/2 - Unp) ^ 4^ x (opp - a„„). 



One can use the above formalism to do consistent and systematic calculations of 
isospin violation. For example, the isospin-violating potential of ranges ~ 



and ~ l/2m7r up to = i^min + 3 were computed in Refs. (|99| , |l00| , |lOl|) . In 
contrast to previous attempts lacking an EFT framework, the EFT results are 
invariant under both gauge transformation and pion-field redefinition and have 
simple forms. The component of range ~ comes from diagrams with all 

possible one-photon dressings of OPE, plus the relevant counterterms (p9|). Its 
isospin structure allows only charged-pion exchange and therefore is class II. This 
TT'j potential has been incorporated in a Nijmegen phase shift reanalysis of np 
data below 350 MeV ([99|). We can use the values for the ttNN coupling constants 
determined by the analysis to find that their isospin breaking (/3i) is consistent 
with zero, with an uncertainty comparable to our expectation from dimensional 



analysis and from n — r] — rj' mixing (102). Similarly, the two contact interactions 



(7s^j) might be viewed as originating in p — u> mixing and pseudovector-meson 
exchange (in particular close- lying doublets such as ai — /i) ( |102| , 103 ). The 



components of range ~ l/2m-,j come from two sources. One is the pion mass 
difference (Am^) in TPE that generates a class II potential (100); the other is 



a TTirNN seagull that arises as a chiral partner of the nucleon mass difference 
(Amjv), and produces a class III TPE potential ( |101| ). All these effects are 
relatively small, with estimated contributions to the scattering lengths of ~ ±0.5 
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fm. These long-range components could be included in the new Nijmegen PSA 

A fit to A^A^ phase shifts including various isospin-breaking interactions was 
carried out in Ref. ( |104| ), improving on an earlier analysis with perturbative 
pions ( |105| ). A slightly different counting is used that enhances by one order 
the electromagnetic effects, electromagnetic corrections being counted as a ~ 
{Q/Mqcd)"^ ■ The EFT potential to order v = Vmin + 2 (including Coulomb) 
is then fitted to low-energy S and P, np and pp phase shifts, using cutoffs in 
the range 300 — 500 MeV. With modest cutoff dependence, isospin breaking in 
the scattering lengths can be accommodated, and higher energies and partial 
waves predicted. A next-order calculation should achieve the level of precision of 
modern phenomenological potentials. 



3.3 Three- and four-nucleon systems 

When few-body systems are considered, one needs to address the issue of few- 
body forces, which, being in general not forbidden by symmetries, will at some 
level contribute to observables. One of the advantages of a field-theoretical frame- 
work is the possibility of the derivation of consistent few-body forces, free of off- 
shell ambiguities. In the standard nuclear-physics approach, few-nucleon forces 
are either inspired by arguments that are independent of the assumptions invoked 
in the NN potential, or simply guessed on phenomenological terms. 

The pionful EFT offers some insight into few-nucleon forces. In addition to 
contact interactions as in the pionless theory, one has further pion-exchange com- 
ponents. The potential, being defined as a sub-amplitude, includes (for A > 2) 
diagrams that have C > 1 separately connected pieces. An n-nucleon force is a 
contribution to the potential that connects n nucleons. 

Weinberg's power counting embodied in Eq. (|4^ ) suggests a hierarchy of few- 
nucleon forces. As in the two-nucleon case, this power counting relies on an 
implicit assumption about the scale appearing in contact interactions. We have 



seen in Sect. 2.2.2 that in the pionless EFT the running of the renormalization 
group toward low energies enhances the size of three-body forces. The latter get 
contaminated by the fine-tunning present in the two-body sector. Whether the 
same happens at the higher energies relevant to the pionful theory is not clear. 
Part of the 3A forces in the pionless theory match onto diagrams of the pionful 
theory which are the iteration of the NN potential (through NN intermediate 
states where at least one nucleon has momentum 0{m-,j)). It is conceivable that 
the enhancement is removed from contact interactions once the pion is introduced 
explicitly in the EFT. 

The new forces that appear in systems with more than two nucleons have been 
derived in Refs. (^,^,^). The dominant potential, at u = 6—3A = i^min, is simply 
the NN potential of lowest order that we have already encountered in Sect. |3.2.3| . 
We can easily verify that, if the delta is kept as an explicit degree of freedom, 
a 3N potential will arise at v = Vmin + 2, a AN potential at v = Vmin + 4, and 
so on. It is (approximate) chiral symmetry therefore that implies that n-nucleon 
forces VnN obey a hierarchy of the type 

{y(n+l)N) Q \ ^ , . 

with {VnN) denoting the contribution per n-plet. (This hierarchy is a non-trivial 
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consequence of chiral symmetry, as there exist non-chiral models that produce 
large three-body forces.) As we discussed in Sect. |3.2.2| , we expect |(1^27v)| ~ 10 
MeV. Using m-,^ and nip for Q and Mqcd respectively, the suppression factor is 
~ 0.05, give or take a factor of 2 or 3. These estimates are in accord with detailed 
few-nucleon phenomenology based on potentials that include small 3A^ and no 
AN forces. This is shown in Table |^ in the case of the AV18/IL2 potential ( [1061 ). 

It proves instructive to look at the form of the first few terms in the few-nucleon 
potential. At = Umin + 2, in addition to corrections to the A^A^ potential, 
one also finds diagrams involving either three nucleons or two pairs of nucleons 



connected via leading contact interactions and static pions. One finds (y,|,10) 
that the various orderings of these diagrams cancel among themselves and against 
contributions from the energy-dependent piece of the iterated NN potential that 
appears at the same order. Alternatively, redefining the potential to eliminate 
energy dependence leads to no 3 A TPE forces of this type at all (|89|). Remaining 
to this order are only 3A forces generated by the delta isobar, if the delta is kept 
explicit in the EFT. At i'mm-|-3 further terms with similar structure arise (11,10) 
— see Fig. (16). 

If, instead, the delta is integrated out, its contributions appear through the 
parameters of the potential at u = Vmin + 3. In this case, there are no 3A forces 
up to u = Umin + 3. In Ref. (107) the u = Vmin + 2 deltaless potential of Ref. 
( |78| ) was used to predict properties of the 3 A and 4A' systems. When the cutoff 
was varied between 540 and 600 MeV, the low-energy 3A^ scattering observables 
showed relatively little cutoff dependence, and came out similar to conventional- 
potential results. The 3A^ and 4A binding energies also fell into conventional 
range, but they changed by about 10%. This nearly- linear dependence on the 
cutoff could be indication of inconsistent renormalization (presumably curable 
by a short-range 3A^ force), although the small range of cutoff variation does not 
allow firm conclusions. It is interesting that when certain v = Vmin + 3 A A 
contributions are included, the cutoff dependence decreases. In the process one 
loses the agreement for the elastic neutron analyzing power, Ay. An "Aj^ puzzle" 
plagues conventional models: potentials that fit AA data well all fall short of 
reproducing Ay at energies as low as 3 MeV, even when existing 3A potentials 
are added. This seems to be the case with chiral potentials as well (108). 



Table 1: Contributions of the two-, three- and four-nucleon potentials (per dou- 
blet, per triplet, and per quadruplet, respectively): Weinberg power counting (W 
pc) and calculations (|106| ) with the AV18/IL2 potential for the ground states of 
various light nuclei (^H, ^H, etc.). 
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^Li 
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~ 10 
~ 0.5 
~ 0.02 
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23 
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0.38 
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? 
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Figure 16: Some time-ordered diagrams contributing to the 3A^ potential in the 
pionful EFT. (Double) solid lines represent nucleons (and/or deltas), dashed lines 
pions, a heavy dot an interaction in C^^\ and a dot within a circle an interaction 
in C^^\ First line corresponds to = Umin + 2, second line to = i^mm + 3, and 
". . ." denote ly > i^mm + 4. All nucleon permutations and orderings with at least 
one pion or delta in intermediate states are included. 





The leading 3A^ potential (10) has components with three different ranges: 
TPE; OPE/short-range; and purely short ranged. (Relativistic corrections ne- 
glected below are discussed in Ref. (0).) The TPE part of the potential is 
determined in terms of vrA^ scattering observables (10). It is similar to the Tucson- 
Melbourne (TM) and Brazil potentials (|109| ), but it corrects a deficiency of the 
TM potential in regard to chiral symmetry. The TM potential was built from 
a vrA^ scattering amplitude that corresponds to a particular choice of pion fields 
for which chiral symmetry is not respected term by term ( |11C| ). Although the 
on-shell TM ttN amplitude agrees with that obtained in EFT to sub-leading 
order, they differ off shell. Of course, in a field theory the corresponding 3A^ 
potentials involve other interactions that enforce chiral symmetry and invariance 
under pion-field redefinitions. Unfortunately the original TM potential was not 
derived within field theory, and the required extra terms were not included. After 
this is done, the EFT result is obtained ( |110| ). The resulting TM' potential is 
similar to the Brazil potential, and has been studied in detail in Refs. ( 112| , 111[ ). 
In Ref. ( |112| ) it was shown that one of the components of the force is dominant in 
3A^ elastic scattering observables. This explains why all existing TPE 3A^ forces 
give essentially the same results for the 3A^ continuum after being fitted to the 
triton energy. This type of 3A^ force does not improve agreement for Ay much 

The novel OPE/short-range components of the potential involve two it{NN)'^ 
interactions of strengths not determined by chiral symmetry alone (^). These 
parameters can be thought of as representing short-range effects such as a and 
uj exchanges with an intermediate A^(1440) resonance, and p exchange from a vrp 
Kroll-Ruderman term ( |113D . It can be shown that due to the antisymmetry of 
the wave- function only one combination of parameters contributes ( |114| ). This 
combination can be determined from reactions involving only two nucleons, as 



discussed in Sect. below. In Ref. (|112| ) the effect of these novel terms was 
estimated assuming that their LECs have natural size. In conjunction with TM', 
this force can bridge a significant part of the discrepancy between "realistic" 
A'^A^ potentials and Ay data. Although not a consistent EFT calculation, this 



Effective Field Tlieory for Few-Nucleon Systems 



43 



estimate serves to assess the size of 3A^ forces in the EFT. Clearly, agreement 
with Ay data at v = in the deltaless EFT is purely accidental. Conversely, 

a full v = Vmin + 3 Calculation might resolve the Ay puzzle. 

The purely short-range components of the potential can only be determined 
from few-nucleon systems (|lO|). As we have seen in Sect. 2.2, the Pauli principle 



leads to a single LEC. Once this LEG is determined from one 3A^ input (such as 
the triton binding energy), all other observables (such as Ay) can be predicted 
(once the OPE/short-range component has been fixed by data involving just two 
nucleons) . 

It is apparent that the pionful EFT brings new forces into play, and that these 
new elements might resolve remaining issues in the description of data. This 
prospect calls for a fully consistent i^mm + 3 calculation with maximal NN input. 

3.4 Processes with external probes 



The power counting arguments of Sect. 3.2.2] can be easily generalized to the 



case where external probes with momenta Q ^ Mjvat interact with few-nucleon 
systems. The probes deposit an energy ~ Mjvtv onto the nuclear system, so that, 
if we define the kernel K as the sum of irreducible diagrams to which the probes 



are attached, the power counting Eq. (43) applies equally well to K. Interactions 



among nucleons occurring before or after scattering can be treated as before: 
iteration of the potential gives rise to the wave- function (|V'')) the initial 
(final) nuclear state. The full scattering amplitude is then 

T = {iP'\K\iP). (59) 

The pionful EFT can also handle scattering at smaller Q, of course, but then Eqs. 
( ^ , |59| ) have to be modified. When the deposited energy is ~ M'^j^/uin — for 
example, when the incoming probes are photons of momenta Q ~ M'^j^/mjsi — 
there can be intermediate few-nucleon states that are reducible, and the break- 



down of T into kernel and wave- functions is more complicated (115). In this 
situation a perturbative treatment of pions, or even better, the pionless EFT, 
should suffice. 

In practice, it is frequently desirable to minimize nuclear wave- function errors 
by using a high-precision phenomenological potential. That this is a good ap- 



proximation is suggested by a comparison (76) between a simplified version of the 
EFT potential of Ref. (|9|) and the Argonne V18 potential (|90|), which show agree- 
ment in most aspects of the wave- function. The cost of this "hybrid approach" 
(P) is the introduction of an uncontrolled error due to a possible mismatch be- 
tween the off-shell extensions of the kernel and the potential. This error can, on 
the other hand, be estimated by the use of several different phenomenological 
potentials of similar quality. 

As with few-nucleon forces, the factor — 2C in Eq. ( |43| ) implies that external 
probes tend to interact predominantly with a single nucleon, simultaneous inter- 
actions with more than one nucleon being suppressed by powers of {Q/MqcdY ■ 
Again, this generic dominance of the impulse approximation is a well-known result 
that arises naturally here. This is of course what allows extraction, to a certain 
accuracy, of one- nucleon parameters from nuclear experiments. A valuable by- 
product of the EFT is to provide a consistent framework for one- and few-nucleon 
dynamics, whereby few-nucleon processes can be used to infer one-nucleon prop- 
erties. More interesting from the purely nuclear-dynamics perspective are, how- 
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ever, those processes where the leading single- nucleon contribution vanishes by 
a particular choice of experimental conditions, for example the threshold region. 
In this case, certain two-nucleon contributions, especially in the relatively large 
deuteron, can become important. Further examination of the structure of the 
chiral Lagrangian reveals that two-nucleon contributions tend to be dominated 
by pion exchange. Indeed, photons and pions couple to four-nucleon operators 
only at 0{Q /Mqcd) relative to pion-exchange diagrams constructed out of the 
leading-order Lagrangian. Thus power counting justifies the "chiral filter hypoth- 
esis" that was put forward to summarize some empirical results on electroweak 
form factors ( |116| ). This "pion dominance" ensures that two-nucleon contribu- 
tions from the EFT in lowest orders tend to be similar to those in phenomeno- 
logical models that include pion-exchange currents. 

Many processes have been analysed in the pionful EFT. Some of those processes 
are extensions to higher energies of the same electroweak processes described in 



Sect. For example, 

• ed — > ed and deuteron form facto rs (|117l 

• ed — > eNN and parity violation (|118| ) 

• np ^ d'y and meson-exchange currents (|IT 

• rip —> d'y and parity violation ( |12C| 

• pp ^ de^Ve and axial currents ( [12 

• p ^He — > ^He e^Ve and solar neutrino production (|12 

• i^d — > INN and solar neutrino detection (|T2 

• fid ^ v^nn and its measured rate ([124| 



• jd ^ ^d and nucleon polarizabilities ( |125 , 115) 

For details, we refer the reader to more extensive reviews ( [iSD and the original 
papers. Here we briefly discuss those processes more germane to the pionful 
theory, involving pions in initial and/or final states. 

3.4.1 vrd^vrd 

This is perhaps the most direct way to check the consistency of EFT in one- and 
few-nucleon systems. For simplicity, consider the region near threshold with the 
delta integrated out. There the lowest-order, z/ = Umin = —2 contributions to the 
kernel vanish because the pion is in an s wave and the target is isoscalar. The 
J/ = lymin + 1 term comes from the (small) isoscalar pion-nucleon seagull, related 
in lowest-order to the pion-nucleon isoscalar amplitude b^^^ . v = Vmin + 2 contri- 
butions come from corrections to ttN scattering and two-nucleon diagrams, which 
involve besides 6*^''-* also the much larger isovector amplitude These various 
contributions to the vrd scattering length have been estimated (§, |126| ). They were 
found in agreement with previous, more phenomenological calculations, which 
have been used to extract 6^''). Partial sets of higher-order corrections have 
been evaluated in Ref. (|127] ) (and in Ref. ( |128| ) for the related, double-charge- 
exchange process). A consistent v = Vmin + 3 calculation of vrd elastic scattering 
is in progress ( |129D . Eventually, a, v = Vmin + 4 calculation might be required 
in order to determine the chiral-symmetry breaking LEG discussed further 

below in connection with lattice QCD. Charge-symmetry-breaking effects were 



considered in Ref. (130). Note that an alternative approach with perturbative 



pions (p. 311) also seems to accommodate the available experimental data. 
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3.4.2 7Wd^7r°(i 

This reaction offers the possibihty to test a prediction arising from a combination 
of two-nucleon contributions and the single- neutron amphtude. The differential 
cross section for a photon of momentum k and longitudinal polarization ei to 
produce a pion of momentum q is, at threshold, 



3k\ f da 



= \Edf + eL\Ld\\ (60) 

q=Q 



where the electric dipole amplitude Eii{k'^) characterizes the transverse response 
and Ld{k'^) the longitudinal response. 

Ed{0) was studied up to = fmm + 3 with the delta integrated out in Ref. 
( |132| ). There are two classes of contributions, according to whether the external 
light particles interact with one or with both nucleons. The one-nucleon part of 
the kernel is given by standard A = 1 ChPT, with due account of P waves and 
Fermi motion inside the deuteron. The neutrality of the outgoing s-wave pion 
ensures that the leading = —2 = i^min terms vanish. The first two-nucleon 
part of the kernel appears at v = Vmin + 2; it comes from a virtual charged pion 
photoproduced on one nucleon which rescatters on the other nucleon with charge 
exchange. These contributions are actually numerically larger than indicated by 
power counting due to the relatively large deuteron size. Smaller two-nucleon 
terms appear at v = Vmin + 3 from corrections in either nucleon. Results for 
E(i{0) up to = i^min + 3 (^32|) are shown in Table ^. They correspond to the 
Argonne V18 potential ( |90[ ) and a cutoff A = 1000 MeV. Other realistic potentials 
and cutoffs from 650 to 1500 MeV give the same result within 5%. The chiral 



potentials of Sect. 3.2.3 are more cumbersome to use, but it has been verified that 
they give results that are similar to other potentials. Two-nucleon contributions 
seem to be converging, although more convincing evidence would come from 
next order, where loops appear. A model-dependent estimate ( p. 33] ) of some 
^ = ^min + 4 terms suggests a 10% or larger error from neglected higher orders 
in the kernel itself. The single-scattering contribution depends on the amplitude 
for jn 7r°n, in such a way that £^^(0) 1.79 - 0.38(2.13 - -EJ+"^) 

in units of W~'^/m^+. Thus, sensitivity to £'0+"^ survives the large two-nucleon 
contribution at v = Vmin + 2. 

We see that working within the EFT yields a testable prediction, Ed{0) = 



(1.8 ± 0.2) • 10~3/m^+ (HI). It is remarkable that for this process EFT gives 



results that are significantly different from tree-level models of the type tradi- 



tionally used in nuclear physics. For example, the models in Ref. (134) predict 



the threshold cross section about twice as large as the EFT. Most of the dif- 



Table 2: Values for £'^(0) in units of 10~^/m-^+ from one-nucleon contributions 
(lA^) up to = z^min + 3, two-nuclcon kernel (2A) at z^ = I'min + 2 and at 
= I'min + 3, and their sum (lA -|- 2A). 

lA 2A lA + 2A 

^ ^ ^min ~l~ 3 V = I'min ~l~ 2 U = I'min ~l~ 3 Z^ ^ l^min ~l~ 3 



0.36 



-1.90 



-0.25 



-1.79 
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Figure 17: Reduced cross section aji = {k/q)a for neutral-pion photoproduction 
as function of the photon energy. Threshold is marked by a dotted line. The 



star at threshold is the EFT prediction ( 132 ), and squares are data points ( [1351) . 
Figure courtesy of U.-G. MeiBner. 



ference comes from one-nucleon loop diagrams: tree- level models tend to differ 
from ChPT mostly by having a smaller Eq^^\ which increases A test 

of this prediction is thus an important check of our understanding of the role 
of QCD at low energies. Such a test was carried out at Saskatoon ( |135| ). The 
experimental results for the pion photoproduction cross section near threshold 
are shown in Fig. ([T7|), together with the EFT prediction at threshold ( |132D . 
Inelastic contributions have been estimated in Refs. ( 135 , 136 ) and are smaller 
than 10% throughout the range of energies shown. At threshold, Ref. (|135| ) finds 
EaiO) = -(1.45 ± 0.09) • 10-3/m^+. While agreement with the EFT to order 

= T^min + 3 is not better than a reasonable estimate of higher-order terms, it is 
clearly superior to tree-level models. This is compelling evidence of chiral loops. 

A further test of the EFT comes from the coherent neutral-pion electroproduc- 
tion on the deuteron. Ed{k'^) and L^ik"^) were predicted to = Vmin + 2 (with 
no new free parameters) in Ref. ( |137| ), for transferred momenta in the range 

— 0.1 GeV^. Because E^^^'^ is not well reproduced at this order, only the /c^ 
dependence can be tested. This reaction was measured at k'^ = —0.1 GeV^ in 
Mainz (138), and values for |£'(j(— 0.1)| and |L^(— 0.1)| were extracted. If they are 
compared with the results from Ref. ( |137 ) simply shifted by a /^^-independent 
amount in order to reproduce the £'d(0) of the u = Vmin + 3 calculation, then 
there is good agreement for |ii^(^(— 0.1)| but |L(^(— 0.1)| fails by a factor 2 ( |138D . 
Because the calculated Lrf(A;^) is not dominated by a single mechanism to the 
extent Ed{k'^) is, it is possible that it suffers from stronger corrections in next 
order. An extension of these calculations to higher order and beyond threshold 
is also highly desirable. 
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3.4.3 NN NNtt 

This reaction has attracted a lot of interest because of the failure of standard 
phenomenological mechanisms in reproducing the small cross section observed 
near threshold. It involves larger momenta of O [y/n^pifiM) , so the relevant ex- 
pansion parameter here is not so small, {jn-,^ j ra^^'^ . This process is therefore 
not a good testing ground for the above ideas. But [jnT^lmfq)^ is still < 1, so 
at least in some formal sense we can perform a low-energy expansion. In Refs. 



( |139| , 140) the chiral expansion was adapted to this reaction and the first few 
contributions estimated. (Note that — contrary to what is stated in Ref. ( |141| ) — 
momenta ~ y^rrvmJv do not necessarily imply a breakdown of the non-relativistic 
expansion, as p^/m^ ~ {m-,^ / mN){p'^ / tun) is still small.) 

Initial attention concentrated on pp ppir^ at threshold. The lowest-order 
terms all vanish, and the formally-leading non- vanishing terms — an impulse term 
and a similar diagram from the delta isobar — are anomalously small and partly 
cancel. The bulk of the cross section must then arise from contributions that 
are relatively unimportant in other processes: isoscalar pion rescattering, TPE, 
and high-order short-range tt{NN)'^ terms. While the first two contributions are 
calculable, the third involves LECs that can only be fitted or modeled. These 
LECs can be thought of as originating from heavier-meson exchange: pair dia- 
grams with a and uj exchange, and a irpuj coupling, among other, smaller terms 
( |142| ). In Ref. (|142| ) it was shown that a large uncertainty comes from the short- 
range features of the wave-function, so a more systematic study has to await 
the development of chiral potentials that are accurate at the relevant energies. 
Other EFT studies of this channel, including attempts to compute TPE, can be 
found in Ref. ( |143D . More problematic is the situation with the threshold cross 
section of other, not so suppressed channels A^A^ dn, — > pmr. In those chan- 
nels the Weinberg- Tomozawa ttttNN term, fixed by chiral symmetry, dominates. 
Wave-function dependence is much smaller, yet a calculation including leading 
and some sub-leading contributions underpredicts the data by a factor of ~ 2 



(144). A calculation including TPE is badly needed. 

Despite these problems, much can be learned from this reaction in the threshold 
region. One example is charge-symmetry breaking. The nucleon mass splitting 
comes from both the quark mass difference and from electromagnetic effects, 
ArriN = 5mN -\- 5mj\f, with 5mj\f = 0{sm^/MQCD) and 5mN = 0{aMQCD /t^)- 
Determining the two LECs 5mN and 5m jv separately is interesting for several 
reasons: coupled to a lattice evaluation of 5m at, it can be used to extract quark 
masses; it can test quark models that evaluate 5m n'-, and it can constrain a pos- 
sible time variation of a because ^He nucleosynthesis is sensitive to Am^y. These 
LECs contribute in combinations other than Am^ to processes involving pions, 
because the two operators that generate the nucleon mass difference have differ- 
ent chiral partners, which involve an even number of pions. Unfortunately, these 
LECs are hard to measure directly in vrA^ scattering. The forward-backward 
asymmetry in np — > dir^, on the other hand, is sensitive to the charge-symmetry 
breaking from these operators, and it has been calculated ( |145| ). Since the asym- 
metry is related to a ratio of amplitudes, some of the uncertainties in the strong- 
interaction physics are reduced. The asymmetry is being measured at TRIUMF 
( |146| ), at a level that could allow for an observation of the quark-mass-difference 
effect. A related experiment, dd air^, which can address the same issues but 
with different theoretical uncertainties, has been proposed at lUCF (|147|). 
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It is possible that some of the problems encountered at threshold stem from 
the smallness of pion s waves, which show poor convergence also in A = 1 ChPT 
— for example, in neutral pion photoproduction on the proton. Indeed, p-wave 
pion production seems better behaved. Ref. ( |14[1| ) calculates the first two or- 
ders of the cross section for the pp — > ppvr" reaction with initial nucleons in the 
spin 5* = 1 state in the direction of the incoming center-of-mass momentum, as 
function of the outgoing pion momentum in the range 0.5 — 1 m,r- With no free 
parameters, good convergence and reasonable agreement with data are found. It 
was also pointed out that other observables would, at that order, be sensitive to 
a combination of 7r(iVA^)^ LECs that affects the leading 3A^ force, discussed in 
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Sect. 3^. In particular, the amplitude for the ^Sq — > (^5*1 — ^ Di)p transition 
which vanishes in leading order, is very sensitive to this LEG. This amplitude 
extracted from the pp pmr'^ data for pion momenta in the range 0.2 — 0.5 m. 
can be fitted quite well with a natural-sized LEG ( |140| ) . This value for the LEG 
can be used in the 3A^ potential to improve the predictive power of the chiral 
potential. 

All calculations of pion-production observables have involved approximations 
necessary to match the kernel and wave- functions. A critical discussions of these 
approximations can be found in Ref. (|148| ) . Issues such as the size of the contri- 
bution of the ttNN cut, not well accounted for in the common approximations, 
need to be better understood. Pion production is clearly wide open for further 
development. 

3.5 Connection with lattice QCD 

The holy grail of nuclear physics has for some time been its derivation from QGD. 
As we have seen, light nuclei are large objects of size 1/Mnn » 1/Mqcd, 
or larger. Dynamics at this scale can be understood within the EFT, and all 
nuclear information is encoded in the parameters of the EFT Lagrangian. These 
parameters, in turn, are fixed by the physics of smaller distances. If the EFT can 
somehow be matched onto QGD at some scale not far below Mqqd, the EFT 
can be used to predict all of traditional nuclear physics. The EFT allows to split 
the quest for the holy grail in two stages, according to the two energy scales. 

At present, the best hope for a solution of QGD in the regime of large coupling 
constant relevant for nuclear physics is by explicit numerical solution on the 
lattice. However, the large size of nuclei make their direct simulation practically 
and intellectually unsound. A more reasonable goal for is to match with the EFT, 
which requires lattices of size not much larger than 1/AIqcd- We are still far 
from this goal, but a few steps have already been taken. 

One obstacle arises from the difficulty in simulating small pion masses. For 



example, Ref. ( 149 ) computes the and ^Si scattering lengths in quenched 
QGD with > 500 MeV. 

The tBtt dependence of nuclear forces comes in explicitly from pion propaga- 
tors in pion exchange, and implicitly from short-range interactions. For illustra- 



tion, in Fig. (18) we exhibit the deuteron binding energy and the ^Si scatter- 
ing length stemming from the leading explicit m^^ dependence in the expansion 
around the chiral limit (|83|). A higher-order two-derivative contact interaction 
was also included and fitted to the triplet effective range. For the physical value 
of the pion mass, one gets the deuteron binding energy to reasonable accuracy, 
Bd = 2.211 MeV (essentially independent of the cutoff R). In the chiral limit 
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Figure 18: The deuteron binding energy (left panel) and the scattering length 
(right panel) as functions of the pion mass that explicitly appears in the OPE 
potential. Implicit pion-mass dependence was not calculated, as other parameters 
were set to their physical values for all m,r- The dots are quenched lattice QCD 
data. From Ref. (||). 



the deuteron is bound by ~ 4.2 MeV. This value is still somewhat small 
compared to f'^/lmj^ ~ 10 MeV, which one might expect to arise in QCD, and 
therefore one would conclude that the deuteron is still weakly bound in the chiral 
limit! This calculation of the explicit dependence agrees with that obtained 
with the AV18 potential (|l50|) with = 0, of B^}^^^^^ ~ 4.1 MeV. The lattice 



id 

data for the triplet scattering length from Ref. ( |14g| ) are also shown. 

While phenomenological models typically can only vary ttItt in OPE, all aspects 
of ttItt dependence can in principle be determined in the pionful EFT. Since the 
pion mass can be varied up to Mqcd, the EFT can be used to extrapolate 
lattice results to realistic values of mj^. It was pointed out in Ref. (|8^ ) that the 
leading (explicit and implicit) 171^^ dependence of nuclear forces can be calculated 

once the chiral-symmetry-breaking LEG is known. That is because the 

leading 771,^ dependence in /^r, qa and m^v is known from ChPT. Unfortunately, 
determination of C^^'^ requires calculation of processes involving external pions 



e.g. vrd scattering, see Sect. p. 4.1 — at high orders and, consequently, precise 



low-energy data. Alternatively, one can imagine fitting c'^^'^ to the lattice data 
themselves. 

Note that in Fig. ( |l8|) we illegally compared the EFT with quenched QCD. 
Most simulations cannot yet be done in QCD itself, but only in quenched or, 
more generally, partially-quenched QCD (PQQCD), where different masses are 
assigned to valence and sea quarks. PQQCD has a different symmetry pattern 
than QCD, a different low-energy dynamics, and thus a different dependence on 
the pion mass. A proper extrapolation of PQQCD data to smaller pion masses 
requires a partially-quenched EFT. Implications of PQEFT to the A^A^ interaction 
are discussed in Ref. (11511). 



4 OUTLOOK 
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4-1 More-nucleon systems 

We have seen that the EFT paradigm has been extensively apphed to systems 
with A = 2 — 4:, for momenta below and above the pion mass. Work remains to 
be done even for those systems, of course. For example, the chiral expansion of 
the pionful EFT still needs to be better understood for A = 2,3; and A = 4 has 
to be studied in the pionless EFT in order to uncover the role of a 4A^ force that 
could appear in LO. 

As these issues get settled, a natural next step for the EFT program is to 
increase A. There have, in fact, been attempts to extend the paradigm to heavier 
nuclei. 

For example, EFT methods are being used to perfect the nuclear shell model 
( |152| ). The goal here is to take some modern potential model and simplify the 
bound-state problem for large nuclei in such a way as to make a numerical solution 
of the Schrodinger equation feasible. This simplification comes about by reducing 
the dimensionality of the original Hilbert space of the shell model, the effect of 
the highly excited states being included into local operators acting on a reduced 
Hilbert space. The reduced problem obtained by "integrating out" the high- 
energy modes can then be solved by standard numerical methods. 

Another approach is to develop an EFT to handle other nuclei that are, like 
the deuteron, particularly shallow (|153| ). Examples are halo nuclei, where the 
separation energy of one or more nucleons is much smaller than the energies 
associated with a core of the remaining nucleons. 

There have also been many EFT- inspired studies of nuclear matter (A — > oo, 
a = 0) and very heavy nuclei. They have all aimed at identifying the relevant 
degrees of freedom and an expansion parameter that can describe physics for 



densities around the saturation density. See, for example, Ref. ( 154 ). 

Whether these approaches prove to be bona fide EFTs (in the sense used in this 
review) or not, the problem remains of deriving the saturation of nuclear matter 
from an EFT adjusted to describe few-nucleon physics. This is a formidable 
problem: the complexity of the necessary resummations of LO operators increases 
rapidly with A, becoming high already at j4 = 5. Lacking the identification of a 
further expansion parameter, we might, as in QCD itself, have to resort to lattice 
simulations. A step in this direction was taken in Ref. (|155| ), where a toy model 
with no- and two-derivative contact two-body interactions was solved (at zero 
and at finite temperature) on a spatial lattice using Monte Carlo techniques, and 
the interaction parameters were fitted to saturation properties. The next step 
involves using EFT interactions determined from few-nucleon systems. 

These attempts, interesting as they are, fall outside the scope of this review, 
and we refer the reader to the original literature. 



4-2 Conclusion 

For the last couple of years, the pionless EFT has been developed and applied 
to two- and three-nucleon systems. Although for the A'^A'^ system in isolation it 
amounts to nothing more than ERT, the full power of the field-theory arsenal 
comes to fruition when more nucleons and/or external probes are considered. 
We have seen that the extension to the 3N system is full of surprises, such as 
the appearance of limit-cycle behavior and of a relevant three-body force. These 
surprises have been turned into successes, and relatively simple calculations yield 
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results of quality not inferior to polished potentials models. Although limited 
in energy, this EFT can achieve high precision for reactions of interest to astro- 
physics, such as — > dj. 

The older pionful EFT is less well understood. There are hints that the ex- 
pansion has finally been identified as an expansion around the chiral limit, but 
higher orders in the expansion must be studied. Among the higher-order terms 
IS a LEG, (7^^""^ that is the main uncertainty in the extrapolation to the chiral 
and heavy-pion limits. External probes might be able to determine this LEG. 
Anyhow, considerable progress has been achieved in the development of the EFT 
A^A^ potential. Isospin-violating effects are a unique virtue of the pionful EFT 
because they are so tightly linked to the pattern of QGD symmetries. There 
remain issues regarding the size of short-range 3A^ forces, but novel longer-range 
3A^ forces naturally appear and can play an important role in nuclear dynamics. 
An assessment of this progress from the perspective of the historical development 
of nuclear potentials can be found in Ref. (|156| ). 

Yet, most nuclei still await us. 
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